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Abstract 



The system of N scalar particles with Grassmann- valued color charges plus 
the color SU(3) Yang-Mills field is reformulated on spacelike hypersurfaces. 
The Dirac observables are found and the physical invariant mass of the sys- 
tem in the Wigner-covariant rest-frame instant form of dynamics (covariant 
Coulomb gauge) is given. From the reduced Hamilton equations we extract 
the second order equations of motion both for the reduced transverse color 
field and the particles. Then, we study this relativistic scalar quark model, 
deduced from the classical QCD Lagrangian and with the color field present, 
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in the N=2 (meson) case. A special form of the requirement of having only 
color singlets, suited for a field-independent quark model, produces a "pseu- 
doclassical asymptotic freedom" and a regularization of the quark self-energy. 
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I. INTRODUCTION 



The Dirac observables for Yang-Mills theory with Grassmann- valued fermions [|T|, 
Abelian and non-Abelian S(2) Higgs models and for the SU(3)xSU(2)xU(l) standard 
model of elementary particles [ffl] and the corresponding physical Hamiltonians have been 
found in a noncovariant way equivalent to a generalized Coulomb gauge, following the scheme 
used by Dirac [E| to do the canonical reduction of the electromagnetic field with charged 
fermions. See Refs. M for reviews of the method and of the program, which aims to get a 
unified description of the standard model and of tetrad gravity in terms of Dirac's observ- 
ables. 

Then, the problem of how to covariantize these results was started. In Ref. @, the 
system of N scalar particles with Grassmann electric charges plus the electromagnetic field 
was described by defining it on arbitrary spacelike hyper surf aces, which give a covariant 3+1 
decomposition of Minkowski spacetime M 4 , following Refs. |8|,[| (see also Ref. , where 
a theoretical study of this problem is done in curved spacetimes). The new configuration 
variables are the points 2 M (r, a) of the spacelike hypersurface S T [the only ones carrying 
Lorentz indices] and a set of Lorentz invariant variables containing a 3- vector ffi(r) for 
each particle [xf(r) = z^(r,ffi(r))] and the electromagnetic gauge potentials A^(r,a) = 
9 qJa^ A*( z ( t > = ( r '^)]> which know implicitely the embedding of E r into M 4 . 

One has to choose the sign of the energy of each particle, because there are not mass-shell 
constraints (like pj — mf pa 0) among the constraints of this formulation, due to the fact that 
one has only 3 degrees of freedom for particle since the intersection of a timelike trajectory 
and of the spacelike hypersurface E T , with Lorentz scalar 'time' parameter r (labelling the 
leaves of the foliation of M 4 with the S T all diffeomorphic to a given E), is determined by 
3 numbers: a = fji{r). 

Besides a Lorentz scalar form of the electromagnetic first class constraints, one has 4 
further first class constraints 7^ M (r, a) pa implying the independence of the description 
from the choice of the spacelike hypersufaces foliating M 4 . Being in special relativity, it is 
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convenient to restrict ourselves to arbitrary spacelike hyperplanes z^(t, a) = x^(r) + 6^(r)cr r . 
Since they are described by only 10 variables [an origin x^(r) and 3 orthogonal spacelike unit 
vectors generating the fixed constant timelike unit normal to the hyperplane], we remain 
only with 10 first class constraints determining the 10 variables conjugate to the hyperplane 
[they are a 4-momentum and the 6 independent degrees of freedom hidden in a spin 
tensor Sjf] in terms of the variables of the particles and of the electromagnetic field. We 
can make the canonical reduction of the electromagnetic field variables to transverse gauge 
potentials and electric fields at the hypersurface as well as at the hyperplane level. 

If we now restrict ourselves to timelike (pi > 0) 4-momenta [the set of particles plus 
electromagnetic field configurations with p% not timelike is of zero measure in the space 
of all configurations], we can restrict the description to the so-called Wigner hyperplanes 
orthogonal to p^ itself. To get this result, we must boost at rest all the variables with Lorentz 
indices by using the standard Wigner boost L^ u (p S) p s ) for timelike Poincare orbits, and then 
add the gauge-fixings 6^(r) — L^ f (p s ,P s ) ~ 0- Since these gauge-fixings depend on p^, the 
final canonical variables, apart p% itself, are of 3 types: i) there is a non-covariant center-of- 
mass variable x M (r) [the classical basis of the Newton- Wigner position operator]; ii) all the 
3-vector variables become Wigner spin 1 3-vectors [boosts in M 4 induce Wigner rotations 
on them]; iii) all the other variables are Lorentz scalars [in the case under consideration 
they are absent after the canonical reduction to the transverse electromagnetic degrees of 
freedom]. Only the 4 first class constraints determining are left. One obtains in this way 
a new kind of instant form of the dynamics (see Ref. |IT1), the Euclidean covariant 1-time 
rest-frame instant form. It is the special relativistic generalization of the nonrelativistic 
separation of the center of mass from the relative motion [H = ^ + H re i\. The role of 
the center of mass is taken by the Wigner hyperplane, identified by the point x^(t) and 
by its normal The 4 first class constraints can be put in the following form: i) the 
vanishing of the total (Wigner spin 1) 3- momentum of the particles plus electromagnetic 
field p[system] « , saying that the Wigner hyperplane Sp^(r) is the intrinsic rest frame 
[instead, p s is left arbitrary, since it reflects the orientation of the Wigner hyperplane with 
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respect to arbitrary reference frames in Minkowski spacetime]; ii) ±yp 2 — M[system] « 0, 
saying that the invariant mass M of the system replaces the nonrelativistic Hamiltonian 
H re i for the relative degrees of freedom, after the addition of the gauge-fixing T s — r w 
[identifying the time parameter r with the Lorentz scalar time of the center of mass in the 
rest frame; M generates the evolution in this time] . When one is able, as in the case of N free 
particles, to find the (Wigner spin 1) 3-vector fj{r) conjugate to p[system}(~ 0), the gauge- 
fixing if ~ eliminates the gauge variables describing the 3-dimensional intrinsic center of 
mass inside the Wigner hyperplane [if ~ forces it to coincide with x%(t) = z^ir, a = if = 0) 
and breaks the translation invariance a \— ► a + a] , so that we remain only with Newtonian- 
like degrees of freedom with rotational covariance: i) a 3-coordinate (not Lorentz covariant) 
%s — \Jr^ s { x s ~ ^t x °) an d its conjugate momentum k s = p s /\f^ f° r the absolute center of 
mass in Minkowski spacetime; ii) a set of relative conjugate pairs of variables with Wigner 
covariance inside the Wigner hyperplane . As noted in Ref. J/J, the noncovariance of the 
center of mass of extended relativistic systems defines a classical intrinsic unit of lenght [the 
M0ller radius p = V ' —W 2 j 'cP 2 = \S\/c\fP 2 determined by the Poincare Casimirs] to be used 
as a ultraviolet cutoff in the spirit of Dirac and Yukawa in future attempts of quantization. 
Let us remark that this ultraviolet cutoff exists also in asymptotically flat general relativity, 
taking into account the asymptotic Poincare charges. 

In this paper we will extend these results to N scalar particles with Grassmann SU(3) 
color charges plus the SU(3) Yang-Mills field. The resulting covariantization provides the 
tools to covariantize the bosonic part of the SU(3)xSU(2)xU(l) model. To complete its 
covariantization, the description of Dirac and chiral fields on spacelike hypersurfaces in 
Minkowski spacetime is needed (this problem is under investigation |]12|| ). 

The final result will be an expression for the physical invariant mass of the system 
scalar quarks plus the color SU(3) Yang-Mills field in terms of gauge invariant quantities 
(covariant generalized Coulomb gauge in the rest-frame instant form of dynamics). This 
pseudoclassical expression is the starting point for defining a relativistic quark model derived 
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from classical QCD. Till now, only the nonrelativistic quantum quark model is available fl3 



and there is no satisfactory derivation of it from QCD, notwithstanding its phenomenological 
relevance. Instead, we have here the full model with relativistic scalar quarks and classical 
color field and the classical reduced Yang-Mills equations with transverse sources analogue 
of the equations OA± = j± of the electromagnetic case in the rest frame [H : if one could 
guess a reasonable solution of the equations for the color field and put it into the invariant 
mass, one would obtain an effective invariant mass for a true relativistic quark model, even 
still with scalar quarks. Since this is not yet possible, we shall limit ourselves to study 
some properties of pseudoclassical mesons (N=2). In particular, we will show that the 
requirement of having only color singlets, realized in a way suited to define a quark model 
without color field, immediately produces a kind of "pseudoclassical asymptotic freedom" 
besides regularizing the quark self-energy. Also some comments on the difficult problem of 
confinement are done. 

In Section II we give some definitions and some results on spacelike hypersurfaces. In 
Section III the Lagrangian of the system on spacelike hypersurfaces is given and the Hamil- 
tonian formalism is developed till the reduction to the Wigner hyperplane. In Section IV the 
Dirac observables of the system are found and the final form of the physical invariant mass is 
given. In Section V the reduced Hamilton equations and then the associated Euler-Lagrange 
equations in the rest frame are given. In Section VI the relativistic quark model in the case 
N=2 (mesons) is defined and it is shown that there is pseudoclassical asymptotic freedom. 

In the final Section there are some conclusions and some comments on confinement and 
the open problems. 
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II. DYNAMICS ON SPACELIKE HYPERSURFACES 



In this Section we will introduce the background material from Ref. J7J needed in the 
description of physical systems on spacelike hypersurfaces, integrating it with the definitions 
of non-Abelian SU(3) Yang-Mills fields 0. 

Let {S T } be a one-parameter family of spacelike hypersurfaces foliating Minkowski space- 
time M 4 and giving a 3+1 decomposition of it. At fixed r, let z M (r, a) be the coordinates 
of the points on £ r in M 4 , {a} a system of coordinates on S r . If a A = (cr T = r;a = {o~ r }) 
[the notation A = (r, f) with f = 1, 2, 3 will be used; note that A = r and A = f = 1, 2, 3 
are Lorentz-scalar indices] and d A = d/da A , one can define the vierbeins 

z^r,a) = d A z^r,a), d s z%- d A z% = 0, (1) 

so that the metric on S r is 

9A6(t, e) = z\(t, ff)T)p,x%(T, a), g TT (r, a) > 0, 

g(r, a) = -det\\g AS (r,a) || = (det\\z^(r, a) ||) 2 , 

j{r,a) = -det \\g n (T,a) ||. (2) 

If 7 rs (r, ct) is the inverse of the 3-metric g?s(j,a) [y ru {r,a)gus{j,a) = Sg], the inverse 
g AB (r, a) of ^(r, a) [g AC (r, a)g^(r, a) = Sj] is given by 



g Tf (r,a) = -{^g T ^](r,a), 

g «( T , a) = Y s (t, a) + [Ig^g^^] ( r , a), (3) 



so that 1 = g (t, a)gg T (r, a) is equivalent to 



7(r, a) 
We have 



</( " T ' SVt(t, ff) - 7"(t, v)grr(T, a)g TS (r, a). (4) 



4(T,3) = {j2-l> t + g 7 rf'zS){T,a), (5) 
V 7 



and 



where 



r]>" = z>i(T,a)g AB (T,a)z» § (T,a) = 

= (l^ + zh fS z»)(r,a), (6) 



l 2 (r,a) = l, l lt (r,&)z?(T,d)=0, (7) 



is the unit (future pointing) normal to S(r) at ^(r, a). 
For the volume element in Minkowski spacetime we have 



d z = z?(t, ajdrdrUp = dr[z^(r, a)/ M (r, B)]J^(t, a)d i o 



'g(T 1 a)dTd- i a. (8) 
Let us remark that according to the geometrical approach of Ref. Jl(|,one can use 



Eq.(|) in the form z£(t,o) = N(r, 5)P(r, a) + N r (r, a)z^(r, a), where N = yjg/'y = 
y/g TT — Y s 9rf9Ts and N r = g T sl ST are the standard lapse and shift functions, so that 
g TT = N 2 + gf§N f N s , g T f = g fS N s ,g^ = N~ 2 ,g^ = -N f /N 2 ,g™ = 7^ + ^, 

The rest frame form of a timelike fourvector p^ is P M = rj^/p 2 ^; 0) = rj^rj^/p 2 , P 2 = p 2 , 

o 

where r\ = signp°. The standard Wigner boost transforming P ^ into p^ is 



ZA,(p,p)=tf(«(p)) 



(Lt , _ ){Pv+Pu) 

Vu "r ^ 2 ~~ 

P P' P +P 2 

1 + u°(p) 



;2 



v = e£(u(p)) = u M (p) = P^/vyP 
v = r e?(u(p)) = {-u r {p)-X - U ^:^ )- (9) 
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u u 

The inverse of L fM u (p,p) is L^ v (p,p), the standard boost to the rest frame, defined by 

L^(p,p) = L v »(p,p) = iA(pi)| pW . (10) 

Therefore, we can define the following vierbeins [the e(f(u(p))'s are also called polarization 
vectors; the indices r, s will be used for A=l,2,3 and o for A=0] 

4(«(p)) = ^(p,p), 

^Hp)) = L A ,(P,p) = V AB V^ b Hp)), 



<£(«(p)) = ty^(«(p)) = 
e o («(p)) = «a(p), 



(11) 



which satisfy 



<(«(p))£(«(p)) = < ? 
<(«(P))4(«(P)) =»7b, 

iT = = «"(p)«"(p) - £ ^(«(p))<(«(p)), 

r=l 

r^AB = e^(w(p))^e^(w(p)), 
Pa^-e^(u(p)) = Pa-^-e^(u(p)) = 0. 

The Wigner rotation corresponding to the Lorentz transformation A is 

( , „ \ 

it^(A,p) = [L(P,p)A- 1 L(Ap,P)]^ = 



1 

i^(A,p) 



R i j (A,p) = (A- 



f :[(A -i,o ((A" 1 )",-!)^- 1 )^ 



(12) 



(13) 



The polarization vectors transform under the Poincare transformations (a, A) in the 
following way 



e?(u(Ap)) = (R- 1 ) r s A^e:(u(p)). (14) 

On the hypersurface E T , we describe the color SU(3) potential and field strength with 
Lorentz-scalar variables A aA (r, a) and F aA ^(r, a) respectively: they contain the embedding 
E(t) — > M 4 and are defined by 

A aA(r, &) = z^(t, a)A a ^z(r, a)), 

F oAb( t , °) = 9aKb{^ v) ~ d B A aA (r, a) + c abc A bA (r, a)A b s(r, a) = 

= z^(t, B)z\(t, a)F aflv (z(r, a)) = z\{t, a)z^(r, a)[d^A av (z(T, a)) - d v A ail (z(r, a)) + 

+Cab c A bfJi (z(T, a))A cu (z(r, a))]. (15) 

We could have written A^ = z^A A = l^Ai + z^A? [zfi are the inverse vierbeins], so to 

get 

A T (r, *) = N(r, a) A^r, a) + N f (r, a)A f (r, a), (16) 
and we could have used Ai(t, a) as the genuine field configuration variable independent 



from the motion of the embedded hypersurface, as suggested in Ref. ||10|| . However, this 
more geometric formulation is equivalent to the simpler one of Ref. J7| [see its Appendix C] 
for spin 1 fields, so that in this paper we shall go on to use A T rather than A\. 

The generators T a of the Lie algebra su(3) of color [A^ = A afl T a ] in the 8-dimensional 
adjoint representation of SU(3) and those T a in the 3-dimensional fundamental one are [c a bc 
are the SU(3) totally antisymmetric structure constants] 

fa = _fa^ (fa )bc = ^ ^« f 6] = c ? 

where the A a 's era the 3x3 Gell-Mann matrices. 
The covariant derivative associated with A afl is 

(D^) ac = 5 ac d^ + c abc A b » = & - A^) ac (18) 
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and the gauge transformations are defined as [F^ = F aflu T a ] 

^ A"(x) = U- 1 (x)A fl (x)U(x) + U-\x)d II U(x) = 
= A^x) + U-\x) (d^U(x) + [A^x), U(x)\), 
F^(x) » pV(x) = U-\x)F llv {x)U{x) = F^(x) + U~\x) [F^x), U{x)\. (19) 

Here U is the realization in the adjoint representation of the SU(3) gauge transformations. 

The scalar particles with Minkowski coordinates x^(t), i=l,..,N, are identified on the 
spacelike hypersurface S T by 3 numbers ff^r), i=l,..,N, by the equation x%(t) = z^(t, rjiij)) 
[so that Xi(r) = z%(t, ^(t)) + zf(T, Vi( T ))vl( T ) an d Vi = ^i( T )]- As shown in Ref. this 
implies that the mass shell constraint pf — m\ = has been solved and a choice of the sign 
of the energy, rji = signp°, has been done. 

In this paper we consider the case of N relativistic scalar particles with the color 
SU(3) charge of each particle described in a pseudoclassical way |L5| [see also Refs. [16 



for pseudoclassicla mechanics] by means of 3 pairs of complex conjugate Grassmann vari- 
ables 9i a (r),9* a (r), a = 1,2,3, which belong to the fundamental representation of SU(3). 
They satisfy 





= o, 


0* 0* I 0* 0* 


= 0, 


n n* _j_ n* n 
u ia"i/3 ~r "ifiUia 


= 0, 



(20) 

and the Grassmann variables of different particles are assumed to commute 

Qiadj/3 = OjfiOia, i 7^ 3 

a* n* n* n* 

a ia a ji3 — a jP°iai 

Qia9*jp = 9jp@ia- (21) 

The color charges of the particles are Q ia (r) = iE^i^WD^^W [Q* a = Qia 
since T at = -T a }. 
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At the quantum level [[L5|] (see also Ref. [|Tj) the Grassmann variables 9* a , 9 ia go 
into Fermi oscillators b\ a , b ia satisfying the anticommutation relations [bi a ,bj /3 ] + = 5ij5 al 3, 
[bia,bj/3]+ = [bl a ,bjp]+ = 0. For each particle there is a 8- dimensional Hilbert space of 
charge states [the charge operator is Qi a = iY^ a /3bi a (T a ) a /3bi/3] with basis | 0j >, b\ a \ 0, >, 
b\ a b\J0i >, &| a &|gfe| 7 | 0j >: the states with k=0, 1,2,3, oscillators transform like a completely 



antisymmetric representation of dimension 



of SU(3). Therefore, the space of charge 



states for each particle transforms like the reducible representation 1 © 3 © 3* © 1 of SU(3) 

, (A 

of dimension J2k=o = 2 3 - 

[k) 

To select the triplet (quark) or antitriplet (antiquark) representation, we shall add to 
the pseudoclassical theory the constraint 



(22) 



As shown in Ref. JL8J, after quantization Ni is replaced by iVj[Aj] = J2a=ob\ a bia — A 



hi — Ai, where hi is the occupation number selecting the 



-dimensional representation 



of SU(3) and Ai is an arbitrary c-number present due to ordering problems. Therefore, 
with the prescription A4 = lor 2, the constraint Ni « becomes the quantum constraint 
iVi [1] I 0i >= or A^[2]|0j >= 0, selecting the triplet or the antitriplet representation 
respectively for particle 'i'. 

With the constraint Ni rs for each i one has QiaQibQic = because it is pro- 
portional to Nf [= in the Dirac strong sense]. Moreover, since one has 52a Qia = 
-e*Jipdl 1 d l5 52a{T a ) a p{T%8 with J2a(T a ) a p(T% s = ±6 a ^ s - \8 a5 5 Pl [valid in the fun- 
damental representation of SU(3)], one gets 

2 



Y,Qla 



■-N? = 0. 

3 1 



(23) 
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III. THE LAGRANGIAN FOR COLORED PARTICLES PLUS YANG-MILLS 

FIELD 



The system of N colored scalar particles plus the SU(3) Yang-Mills field is described by 
the action 



S = J drd 3 aC(T, a) = J drL(r), 
L(t) = J d 3 aC{T,a), 



■ N 3 3 

C(t, a) = - ffi(r)) Ei d Ur)0ia(r) - C(r)^(r) + \(r) £ 9* a (r)e ia (r)} - 

i=l a=l a=l 

N I : 

- S - Vi^MVi^iy 9rr(r, a) + 2# rf (r, a)fil(r) + g H (r, £)?£(t)7)?(t) - 

i=l 

- Qia{A aT {r, a) + A af {r, a)T)f(r))] - 

a 

-L y^J)g A0 (r, a)g^(r, a) £ F aAS (r, a)F aCb {r, a), (24) 



where the configuration variables are z^(r, a) A aA (r, a), ^(r), 9 ia (T) and 0* a (r), i=l,..,N. 
The particles have Grassmann-valued charges Q a i(T~) = i X^ =1 0* a (r)(T a ) Q/ j0j0(T). 
We have 

-\\J~99 AC 1 g BD T, a F a AB F aCb = 
= -\V9^ a [2{g TT g f ~ s - g Tf -g TS )F aTf .F arS + Ag f ~ s g™F aTf F a - sil + g™g' s *F af - s F a ^ = 

= — /^Y [I Pip . 7 «P . _ fin _^vfp _ sup , _J_ I [In/fsp _p / UV , 



oln -ry^n -V" 3 ")! = — ^-(F ■ — F~ A~/ fS ( F - — N* F~) — N ^ F-- F- 

^ gijTva J iSrn } }\ 2N \ rr ur ) I \ TS vs J 4 ; / ± ru*- sv 

The action is invariant under separate r- and a-reparametrizations, since A aT (r, a) trans- 
forms as a r-derivative; moreover, it is invariant under the odd phase transformations 

The canonical momenta are [E af = F afT and B af = \tfsiFasi ( e rst — eVst ) are the electric 
and magnetic fields respectively; for g^ B ~^ VAb one S e ^ s = ~F a f = E^} 

y/grrir, 3) + 2g T r(r, a)-qt{r) + g fi (T, ^f(r)^00 
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+^^\{9 TT z T , + g Tf z f ,)(r, a)g AG (r, a)g 6£, (r, a) £ F aAS (r, a)F aC£) (r, a) - 
- 2[z t ,(t, a)i !l A: g :r g" n + // U V ; // ")(;. a) + 

+ z f ,(r, a)(g Af g^ + g A ^g fd )(r, a)g Si) (r, a)] £ F aA6 (r, B)F a66 (r, a)] = 

a 



dd T A aT {T,a) 

<(t, B) = dL = T^2L 7 «( r , a)(F arS + g T ^F a , s )(r, a) = 

= 9s - 2 ^(r, B)(E aS (r, ff) + g T v(r, a) 7 M (r, ^c^S^r, a)), 

vM r > CT ) 



Kifij) — 



9L(r) _ 
drjl(r) 

g Tf (T, ffi(T)) + g fi (r, ^(r))?)|(r) 



TjiUli 



\jg~^{T r Mj)) + 2^ T r(r,7/i(r))?)[(r) + 0„(r, ijiM^SW 
- I]Qia(r)A^(r,?/i(r)), 

a 

^^^^^W' (25) 
and the following Poisson brackets are assumed 

{z^r,a),p u (r,a'} = -^5 3 (a-a'), 
{A aA (r, a), Trf (r, a')} = r,f5 ab S 3 (a - a'), 

{0ia(r),ireji3(T)} = —8ij8 a pi 

{0Ur),^ 3P {r)} = -5 l3 5 aP . (26) 
The Grassmann momenta give rise to the second class constraints n dia + ps 0, 
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vi> ia + \0i a ~ [{-K eia + |0* a , TT e * j/3 + \Q jfi \ = -iS i:j 8 a p}; TT eia and ir e * ia are then eliminated 
with the help of Dirac brackets 

{A, BY = {A,B}-tj2 [{A, 7r eia + l -e* a }{7re* ia + \e ia , B} + 

a=l 

+ {A,7T0* ia + % -e ia }{Keia + ^0* a ,B}] (27) 

so that the remaining Grassmann variables have the fundamental Dirac brackets [which we 
will still denote {., .} for the sake of simplicity] 

{0 ta (r),9 jP (r)} = {9Ur),9*,(r)} = 0, 

{6 za (T),0*p(T)} = -t5 tJ 5 a p. (28) 

These equations imply {Qi a , Qjb} = 8ijC a bcQic and that the gauge transformations S9i a = 
(a a T a ) a p9ip, under which the action is invariant, are generated by the QiaS. 

By varying the Lagrange multipliers Aj(r), we get the Grassmann constraints mentioned 
in Section II 

^) = ECW9«W»0. (29) 

a=l 

We could also treat the multipliers Aj(r) as configuration variables: we would get the first 
class constraints ^(r) ss and we could get free of the Aj's by adding the gauge-fixings 
Aj(r) f« and going to the Dirac brackets for the resulting 2N second class constraints. 
/,From the expression of the momenta we obtain the four primary constraints 

T-tix(r, B) = p M (r, a) - / M (r, a)[T TT (r, a) + 

N 

i=i 

Vi J m i -Y' s {T,a)[Kif{r) + ^Qia{r)A a f{T,a)}[Kis(T) + ^ Q ib (T)A bs (T, a)} ] - 

V a b 

N 



'r,a)Y S (r,B){-T TS (r,a) + £ <5 3 (<7 - ^(r))^ + £ Q^K^t, <?)]} » 0, (30) 



where 
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T TT (r, B) = -\ Y,{^>K f a g»< ~ ^T^Kr, 5% 

* a 9 s 

T tS {t, a) = - F aSi (T, a)vr* (r, B) = -e s{ii ^ tt*(t, B)B aii (r, a) = 

a a 

= J2[7? a (T,a)xB a (T,a)} s , (31) 

a 

are the energy density and the Poynting vector respectively We use the notation J2 a (^a x 
B a ) s = (E a x B a ) § because it is consistent with egiuYla^aBau m the flat metric limit g^g — > 
VAS ; in this limit T TT - \ £ a (£ a 2 + B\). 

Since the canonical Hamiltonian is (we assume boundary conditions for the electromag- 
netic potential such that all the surface terms can be neglected; see Ref. [[L|) 

N , 

H c = -Y. k»(t)%(t) + / d 3 *[J2 K{t, B)d T A aA (r, B) - p M (r, B)z»{r, B) - C(r, B)) = 

8=1 J a 

= [ d 3 aJ2[d f (K(r,a)A ar (r,a)) - A ar (r, a)T a (r, a)] = 

= - [d 3 a^2A aT (r,a)T a (r,a), (32) 
with [d f = £, = -d\ d = {&}, A - d 2 , D ab = {5 ab d f + c acb A f c } ] 

N 

r o (r, B) = <9r<(r, a) + c abc A bf (T, B)tc t c (t, B) + Qia(r)5 3 (a - 7fc(r)) = 

i=i 

\(-A) X 



= -D a6 (r,(r)-7r fe (r,(T) + ^g w (r)5 3 (a-^(r)), (33) 

i=l 

we have the Dirac Hamiltonian (A M (r, <?), A r (r, a) and Hiij) are Dirac's multipliers) 
# D = J d 3 a[X^(r, B)H,{t, B) + Y. A -(r, a)<(r, a) - £ A aT (r, a)r a (r, B)] + 



iV 



+ ^^(r)iV i (T). (34) 
i=i 

The Lorentz scalar constraints vr^(r, a) ~ are generated by the electromagnetic gauge 
invariance of the action; their time constancy will produce the only secondary constraints 
(Gauss laws) 

r o (r,*)«0. (35) 
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The six constraints Ti^ij, a) ~ 0, vr^(r, a) ~ 0, T a (r, a) « are first class with the only 
non vanishing Poisson brackets 



Moreover, since {Qi a , Nj} = 0, also the constraints iVj ~ are first class and, being constants 
of the motion for each i, they do not generate secondaries. 

Let us remark that the simplicity of Eqs. (p6|) is due to the use of Cartesian coordinates: 
had we used the constraints Hi(r,a) = P(r, ^TC^t, a), Hf(r,(r) = Zf(T,a)H^(r,a) (i.e. 
nonholonomic coordinates), so that their associated Dirac multipliers Xi(t, a), A^(r, a) would 
have been the lapse and shift functions of general relativity, one would have obtained the 
universal algebra of Ref . || . 

The ten conserved Poincare generators are 



so that the total momentum is built starting from the existing energy momentum densities 
on the hypersurface 



{n^r,a) , H u (t,ct')} 





(37) 




N 



+ £5 3 0?-^(t))x 



i=l 



rji mf -j fS (r, a)[K if {T) + ^Qi a (r)A a f(r, a)][K is (r) + ^2Qib(r)A bs (T, <?)]] 





N 



+ £ 6%a - Ur))Ws{r) + £ Q ta (r)A aS (r, a)}} « 0. 



(38) 



i=l a 



We add the gauge-fixings to arbitrary hyperplanes as in Ref. |?J 
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C(t, a) = z»(t, a) - a£(r) - b^r)a f » 0, (39) 

where bf(r), f — 1,2, 3, are three orthonormal vectors such that the constant (future point- 
ing) normal to the hyperplane is 

P( r , B) » P = 6>? = e^ 7 &?(T)&?(r)&2(r). (40) 

Therefore, we get 

;#(t,*)«#(t), 
;#(t,^«#(t) + #0V, 

^«(r, 5) « -£«, 7 H (t> 5) « -<f s ", 7(r, a) « 1. (41) 

By introducing the Dirac brackets for the resulting second class constraints [now we have 
K(r),«5(r)} = ( 5y5«] 

{A, B}* = {A, B}-J d 3 a[{A, C(r, a), 5} - {A, Tt^r, a)}{^(r, a), B}], (42) 

we find {^(r), ^(r)}* = —77^ [with the assumumption {bf{r),p v s } = 0]. 

The ten degrees of freedom describing the hyperplane are x%(t) with conjugate momen- 
tum and six variables 4>\{j), A = 1, .., 6, which parametrize the orthonormal tetrad b^(r), 
with their conjugate momenta T A (r). 

The preservation of the gauge- fixings C M ( r ) (?) ~ in time implies 

^rC M (r, (?) = {CXr, (?), H D } = -A"(r, (?) - a* (r) - ftjfOV » 0, (43) 

so that one has [using = and 6£(t)6£(t) = -&£(t)&{?(t)] 

A^(r )( ?)^^(r) + A^(r)6^(r)^, 
A>(r) = -x^r), 

>T{r) = -A""(r) = ^(r)^(r) - 6j?(r)#(r)]. (44) 
Thus, the Dirac Hamiltonian becomes 
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1_ > N 

H D = A"(r)«,(r) - -A""!^,^) + ^X^A^r), (45) 

i=l 

and only the following first class constraints are left [now we remain with the variables 

h% S^, A A , 7T A , ff h K h 9 ia , 9* a ] 



J d 3 aH^r, a) = jf s - Z„{± J dV^[^(r, a) + g; 2 B 2 a (r, a)} + 



N 



+ E ViJrn 2 + [k 2 (t) + Q ia (r)A a (r, jfc(r))] } - 

i=i V a 

- 6£(r){ / rfV^ [7r a (r,a) x 5 a (r,<7)], + EM^) + E Qia(r)A a ,(r, jft(r))]} « 0, 
7#"( T ) = #(r) / dWW(<r,*) -&£(r) J d 3 aa f W{r,a) = 

= S?(t) - [#(r)6£ - Vi{T)b»\ [i / dW El9>1(^ <?) + ^ 2 5 a 2 (r, a)] + 



iV 



+ E^kW + ftM +E^(^a(r,^(r))] 2 ] + 

i=i V a 

+ mr)Ui{r) - 6j;(T)6ff(r)] [ / dW E [#a(r, a) x B a (r, a)] s + 

N 

+ E^ f WK : «+EOi«(^(^.W)] ] « °' 

i=l a 

<(r, a) w 0, 

r a (r, a) «0, 

iV, « 0. (46) 
Here S^" is the spin part of the Lorentz generators 

= #(r) J d 3 cra f p u (r, a) - b v f {r) J d 3 aa f ff(r, a). (47) 

As shown in Ref. ]7| instead of finding A (r),T A (r), one can use the redundant variables 
6^(r), S'^(r), with the following Dirac brackets assuring the validity of the orthonormality 
condition if - b^rj M b u B = [C$ aP = VjVsV^ + VftsV™ ~ VyVsV^ ~ VyVsV up are the 
structure constants of the Lorentz group] 
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{sfry A y = rrvA-rrvA 



The Dirac brackets of the left constraints are 

{W"(r), W"(r)}* = £ / d 3 a{m(r) ~ b^(r)]g 2 s n af (r, a) - 
-b^r)F ars (r,a)bUr)}T a (r,a), 

{W"(r),W^(r)r = -S; / rfVa*{[(6 r "(r)6f(r)-^(r)^(r))6^ - 

a J 

-(b^(r)-b^ T bUrm(rM^(r,a) + 

+6^(r)F a , l (r,o : )[&f(r)6f(r)-6f(r)6?(r)]}r a (r,a), 
{7i^(r),7i; a/3 (r)}* = C^H< s {t) + £ / rfVaV* 
{6S(r)^(r)([6?6?(r) - ^(t^tt^t, a) - b v f {r)F af - s {r, a)bl{r)) - 
-6g(r)6j(r)([^(r) - «t)]<? s 2 tw(t, a) - ^(t)F^(t, (r)) - 
-bl(rmr)([b^(r) - b^(r)]g 2 s ir af (T, a) - ^(r)F^(r, <?)&? (r)) + 
+6£(r)6j(r)([^(r) - b a T b^r)]g 2 s n af (r, a) - ^(r)i^(r, *)&«(r))}r a (r, a), 



(48) 



(49) 



Let us now restrict ourselves to configurations with p s > and let us use the Wigner 
boost L^ u (p s ,p s ) to boost to rest the variables b^, Sg" of the following non-Darboux basis 



x»,p»,b%Sir,rti,4 



of the Dirac brackets. The following new non-Darboux basis is obtained [x% is no more a 
4-vector] 

~n a , 1 A ( ( de%(u(p 8 )) 

x s = x s + 2^W«))^s— ^ 5f = 



i 



+ 



+ 2i] s Jp 2 s i] 



P°s+VsJp 2 s VsJp 2 s(P°s+VsJp 



■S?p r s }, 
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hi = <(«(p.))#, 



sr 



jr = ~<p v s - viz + sir = v? + sr. 



We have 



{5?,rfr = -ir, 



Ps + Vsypl 

{Si j ,b r A }* = (S ir S js - 8 is 8 jr )b% 

J Q^ v QCtf}\* Q7<5 

\°s ' s J °7(5 °s J 

As shown in Ref . J7| , under Poincare transformations (a, A) we get 
p'* 4 = A"/, 



A^K + - 

A°aP? + V^P 2 s 
(Ps + VoV\ffs)Ps,rA 



(A° Q - l) Ps 

p° s + vJpI 



V\/P 2 s(P°s+V\/P 2 s) 

Therefore, is not a 4-vector: its infinitesimal transformation properties under 
transformations generated by J£ v are 
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{a* J?} = rTSZ ~ V^^ + Sfh 



p° + vyp 2 vyp 2 s(p°s + vyp 2 s ) 

= 0. (53) 
We can define [the new variables are x^p^ 6^, S% v , A aA , tt^, ff i: K i: 9 ia: 9* a ] 

« [tf W " bf(r)bf] [\ J d 3 aa f &) + 97 2 B 2 a (r, *)] + 
N I ~ 

+ E^W^K + [Zi(r) + y EQ i a{T)Mr i fk(r))] 2 ] - 
i=i V a 

- [bHr)bf(r) - bf{r)bj{r)\ [ / E [*a(r, a) x B a (r, a)] s + 

N 

+ E<(r)[<(r) +EOia(r)^(r,^(r))] ]. (54) 

i=l a 

Let us now add six more gauge-fixings by selecting the special family of spacelike hy- 
perplanes orthogonal to p^ (this is possible for p 2 s > 0), which can be called the 'Wigner 
foliation' of Minkowski spacetime. This can be done by requiring (only six conditions are 
independent) 

t a(t) = Ufa) - <U(«(p.)) « o 

=► bi(r) = e>(p s ))^(r) « rjj. (55) 

Now the tetrad has become e^(u(p s )) and the indices £ f ' are forced to coincide with the 
Wigner spin-1 indices 'r', while o = r is a Lorentz-scalar index. The final Wigner- covariant 
variables are x%, p%, A ar , -n T a , A ar , n r a , rjl, n ir , ia , 6* a . One has 

Sf B - (Vfv? - V?t£) [\ J <*W U9st(r, a) + g ; 2 B 2 a (r, a)} + 



N 
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i=l 



J a 

N 

+ E^)[«!W+E«»( T K( T .M ]» 



i=i 
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N 

i=l a 
AT 

- E^( T )K( r ) +EQk(rK(T,#))] ) + 

i=l a 

+ / rf3(T [^a(r, a) x 5 a (r, a)]" - a s [7r a (r, a) x S a (r, a)f ) 



A' 



§f « -£f = -£i£(t)i* /m? + [Kj(r) + EQ»(r)A,(r,ifc(r))]2 

i=l V a 



.r « ,-// - 1*±=. (56) 

P° s +Vs\Pj 



The time constancy of Tjf ~ with respect to the Dirac Hamiltonian gives 
4-mr) - ^(«(P.))] = {tf(r) - ffj,}' = 



l -\ a ^r)mr),S sa p{r)Y = A^(r)6 ?a (r) « 



A^(r) ss 0, (57) 

so that the independent gauge- fixings contained in Eqs.(|55"D and the constraints IrV^ir) « 
form six pairs of second class constraints. 

Now we have [remember that (r) = — A M (r)] 

F = K = u"(p a ), 

z»{t) = x?(t) = ^(p,)-M#(P.)K(«(P.)). 
»(t) = [^(r)^(p s ))] 2 , 

# rr = ^ rr = -, g Tf = -x sll 5 fs e»(u( Ps )), 

9 9 

9rf = is^ s e^(u(p s )), g rs = -5 rs + x slx e^{u{p s ))x sv e v v {u{jp s )). (58) 

On the hyperplane £w r all the degrees of freedom z M (r, a) are reduced to the four degrees 
of freedom x%(t), which replace x^. The Dirac Hamiltonian is now Hp = A M (r)7^ A1 (r) + 
Eilx ^(r)N t (r) + J d 3 a £ a [A ar (r, <?K(t, a) + A a (r, <x)r a (r, a)] 
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To find the new Dirac brackets, one needs to evaluate the matrix of the old Dirac brackets 
of the second class constraints (without extracting the independent ones) 



c 



{H afS , W & }* « {H a(3 , T|}* = 

= S 6B [r,^(u(p a )) - V aa 4WPs))} 
{T P A ,W 5 Y = {Tj,TJ}*=0 

\ = 5 AA \ri pl t S AHPs))-ri p5 t\nPs))} 



(59) 



Since the constraints are redundant, this matrix has the following left and right null 
/_ \ / ~ \ 



eigenvectors: 



a a (3 — ap a 




[a a/ 3 arbitrary], 







£{u(Ps)) 



. Therefore, one has to find a left 



and right quasi-inverse C, CC = CC = D, such that C and D have the same left and right 
null eigenvectors. One finds 

/ 

C = 



\ h,r ef ( Ufa) ) - t] 5t e°(u(p s ))\ 

{ zKv^Hps)) - v^vHps))] o* t d 



CC = CC = D = 



V 



QocfW 



Urfiril - e D "(u(p s ))e AT (u( Ps )) 



(60) 



and the new Dirac brackets are 



{A,B}** = {A, B}* - -[{A^Y^efHPs)) - (u(p s ))]{T^ B}* + 



+ {A,T°Ylv*^Hp s )) - Va ,e«(u(p s ))]{H^,BY] 



(61) 



While the check of {7i a/3 , BY* = is immediate, we must use the relation b A/J Tj^e Dp = —T? 
[at this level we have T\ = T%] to check {T P A , BY* = 0. 

Then, we find the following brackets for the remaining variables x^ : p p , r\\, k\ [the metric 
y s = —§ rs will be used, so that {r/[, .} = d/dn^ = —d/dK ir ] 



and the following form of the Poincare generators [{L^ v , S pv Y* 7^ 0] 



(62) 
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77^ _ ^•t 1 r i 1 ' _ _|_ CM" — T ^ u a_ 

J s — x sPs x sPs +-J S — JJ S "I" O s , 

rir qrs n s 
noi _ s f*g 

§ij = girgjsgrs (63) 

Therefore, is not a fourvector and jfj, «i transform as Wigner spin-1 3- vectors. Indeed, as 
shown in Ref. J7|, under global Poincare and infinitesimal Lorentz transformations one has 

Vir = rjisR'riKp), 
K ir = Ki s R r (A, p), 

$ iS {p r sVi ~P S sVi) 



for, jy'} = 5 is ^'*{ry[, sf } = (<W r - s ir s i ')Tj?, 

r r roi X _ ^'(pX ~PX) 
l K i ) J s J — r— ■ ) 



{4, jv} = (S is 8 jr - 8 ir 8 js )Kl. (64) 
The only left first class constraints are 

W"(r) = - u"(«(p.)) [1 / dVEK(r, a) + g; 2 Bl(r, a)] + 



iV 



+ E W m ? + [^( r ) + E <5ia(r)A a (r, ^(r))] : 
i=i V a 

- e?(u(p 8 ))[[ dVE [vfa(r,a) x 5 a (r,a)] r + 



A/ 



<(r,a) «0, 



r«(r,*)«0, 
iVi(r) « 0, 



{^,^}** = E / <M[eS(u(p.)K(u(p.)) - e^( W (p s ))e« W (p s ))]7r ar (r, 3) + 
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+e?(u(p s ))F ars (T, ay s (u( Ps ))}T a (T, B) 



(65) 



or 



N 



n(r) = - mj + fair) +J2Qia(r)A a (r,ff i (rW + 

i=l V a 

+ \j ^EfoMfa *) + g; 2 B 2 a (r, b)\ } » o, 

JV 

8=1 a 

+ J d 3 crJ2 ^(r, a) x S a (r, a) » 0, 



<(r,a) w 0, 



r a (r,a)^0, 
^(r) w 0. 



(66) 



The first one gives the mass spectrum of the isolated system, while the other three say 
that the total (Wigner spin 1) 3- momentum of the N particles on the hyperplane 
vanishes. The Dirac Hamiltonian is now Hp = X(t)TC(t) — A(r) • TC p (t) + Y,f=i A i i( r )^Vi(' r ) + 
/ d 3 a£ a [A aT (r, a)<(r, B) + A a (r, B)T a (r, B)\ and we have £J = H D }** = -X(r)u^{p s ). 
Therefore, while the old had a velocity x^ not parallel to the normal / M = u^(p s ) to the 
hyperplane as shown by Eqs.(|58|), the new x{f has and no classical zitterbewegung. 

Moreover, we have that TL = I ■ aL = I ■ x„ is the Lorentz-invariant rest frame time. 



For N > 2, let us perform the following canonical transformation [a 



N-l] 



X? Ps 







T 

± s 


z s 


k s 










-1 






Pa 






V+ 







(67) 



Ps ' %s Ps ' %s 



Vs\/P 2 s Vs\/P 2 s 
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whose inverse is 



N 

K + = 51 K i j 
i=l 
_ N 

Pa = ^/N^JaiVi, 
i=l 

1 ^ 



a VNtl 



N N 

i=l i=l 
JV-1 j 

^ lailaj — $ij T7 ; 
a=l iV 



x°. = y/T+%(T. + —), 
x s = h (T s H )/c s , 

Ps ^sksi 



I JV-1 

Vi = V+ + ~7= £ 7ai/?a> 
ViV s=1 

1 JV_1 
iV a=l 



The new form of the constraints is 
n(r) = e. - Q2th x 

i=l 



j JV-1 JV-1 

m 2 + [ £ + ( r ) + v/iV E 7«v?a(r) + £ Qia(r)A a (r, ff + (r) + -= ]T 7«Pa(r))] 2 

\ 7V a=l a VN g=1 
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+ \j d 3 aJ2[9 2 X(r, 3) + 97 2 B 2 a (r, a)}} = e s - E (P+I)s - E (F)s » 0, 



N ^ N-l 



Hp{T) = K+(t) +J2Y,Q™( r ) A a( T >V+(T) + -j= J2 %iPa( T )) + 



i=l a VJV a=l 

+ M 3ff E ^( T > x ^«( r ' <?) = + ~ o, 

J a 

<(r,a) « 0, 
r a (r, a)«0, 

7V.( r ) » 0, (70) 

where £ (F)s = \ J d 3 a^2 a [g 2 nl{T, a) + g~ 2 B 2 a {r, a)] and P {F)s = J d 3 a £ a ^(t, a) x B a (r, a) 
are the rest-frame field energy and three- momentum respectively [now we have n a {r, a) = 
g~ 2 E a (T, a)], while E( P+I ) S and P(p+i) s denote the particle+interaction total rest-frame en- 
ergy and three-momentum, before the decoupling from the electromagnetic gauge degrees 
of freedom. 

The final form of the rest-frame spin tensor is 

N i N-l s ( T \ N-l 

i=i viV s= i iv 5=l 



+ E^K(t,77 + (t) + -= E «r))) - 

a V iV c=l 



V^V s=i ^ i=1 



N-l 



+ J2Qia(r)A r a (r : ff + (r) + -= E 7«Pc(r)))} + 

/ d 3 a EK [*«(r, 5) x B a (r, a)}* - a s [* (r, a) x fl ( T , a)f) = S£ +J> + Sfo, 



+ 



nor nro 

D s — — D s 



V 1 N-l 



= -EWW + 7f E TWSOOH x 
i=i V iV a=1 



i JV-1 i N-l 

\ ™ 2 + [-K + (r) + v / iVE 7»v? a (r) + E Qia(r)A a (r, ff + (r) + -= E 7^(r))] 2 
N 7V a=i a VN - a=1 

\ J d 3 aa r U&Kt, a) + g; 2 B 2 a (r, a)], (71) 
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while the Dirac Hamiltonian is 

H d = \(t)H-X(t)H p + f d 3 aJ2[Kr(r,ff)K(r,ff)-A aT (r,a)T a (r,a)] + 

N 

+ $>(r)iV l (T). (72) 
i=i 

On an arbitrary spacelike hypersurface or on the Wigner hyperplane one has in the free 
case [but also in the interacting one looking at the four constraints 7Y M (r) rs of Eq.fl65|) 
and of Ref. 0] 

Pikr( T ) = Viy™[-YHr, a)K is (r)K is (T)l^(T, a) - z£(r, v)Y s (r, o)ku{t) 

Pilwignerhyperplaneir) = Tji^j Vflf + ft? (V) W^(> s ) + e£(u(p,))K?(r), 

which are solutions of — = 0. 
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IV. THE DIRAC OBSERVABLES 



In this Section we will make the canonical reduction with respect to the SU(3) Yang-Mills 
gauge transformations. From now on we shall use the notation {., .} for the Dirac brackets 
{., .}** on the Wigner hyperplane £^ r . 

The decompositions in the non-Abelian SU(3) gauge potential can be obtained from the 
second paper in Ref. M. For the vector potential we use Eqs.(4-13), (4-16), (4-26), (4-29), 
(4-30), (4-31), (4-33), (5-21), (5-24), of that paper to get 

A a (r,a) = A ab (^ A \r,a))d V i A \r,a) + (P e^^^UA^r, a) = 
= e a ( V ^(T,a),dr]^(T,a)) + (P e n ^^ A) ^) ab A b± (r, a), 
d-A a± (r, ff) = Q 



f a A ab ( V ^(r, a))d4 A \r, a) ■ da = H b (^ A \r, B))B V {A \t, a) ■ da = 
= e a {^ A \r,a),d^ A \r,a))f a = ^^(^(r, <?)), a = a • da, 



^\v iA) (r,a)) = Q ( £ > (v iA \r,a))f a 

rrf A \T,a,s) 



(7) 



f T ' CT,S H b ( V W(T,a;s))T>4 A \T,a;s). (73) 
Jo 



If T) a are coordinates in a chart of the group manifold of SU(3), the matrices A ab (t]) 
satisfy the Maurer-Cartan equations, which can be written in the zero curvature form 
9H q^ — 9I q ^ + lHa(v)i Hb(v)] = 0- We shall use only canonical coordinates of the first 
kind, defined by A ab (r])r] b = rj a [ so that A{rj) = with (Tr]) ab = (T c ) ab i] c = c abc r) c ). 

If 9 a = A ab {jf)dri b are the left-invariant (or Maurer-Cartan) one-forms on SU(3), the ab- 
stract Maurer-Cartan equations are d6 a = —\c abc 6 b A 9 C ; then, by using the preferred line 
7^(5) (s is the parameter along it) defining the canonical coordinates of the first kind in a 
neighbourhood of the identity I of SU(3), one can define d^u^ (r)(s)) = 9 a (r](s)) [d( 7ri ) is 
the exterior derivative along j v ] with ^^(^(s)) = (r)(s))T a = ( 7 ) T a A ab (fj)df] b = 
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(t»i) Jj eta ^ 9 a \ lr T a = ( 7 ) Jj v ^ ^g\^i where uj g = 9 a T a is the canonical one-form on SU(3) in 
the adjoint representation. In our case of a trivial principal SU(3)-bundle P(E(r), SU(3)) 
over the spacelike hypersurface S T [diffeomorphic to R 3 } of Minkowski spacetime, it is shown 
in Ref. Q that Q a (r)( A '(T, &) , drj^ (t, a)) and Q^(t]( a >(t, a)) are just the extension of these 
SU(3) objects: in the second paper of Ref. PJ, a connection-dependent coordinatization 
(r, a;r]^ A \r, a)) of the principal bundle is given with the SU(3) fibers parametrized with 
parallely transported (with respect to the given connection) canonical coordinates of the 
first kind from a reference fiber over an arbitrarily chosen origin in R 3 . The functions 
Va ( r > an d their gradients dr)^' (r, a) vanish on the identity cross section 07 of the trivial 
principal bundle. The path 7 is a surface (in the total bundle space) of preferred paths, 
associated with these generalized canonical coordinates of the first kind, starting from the 
identity cross section 07 till a cross section parametrized by the parameter s, in a tubolar 
neighbourhood of 07. The operator d(*j\ is the exterior derivative on the principal SU(3)- 
bundle total space restricted to 7; it can be identified with the vertical derivative on the 
principal bundle and with the Hamiltonian BRST operator. With these conventions, one has 
{.,r Q (r,0)} = {.,-3 ab ^ b (r,a)} = -B ba (v (A) (r,a))-jJ— [B{rj) = A'^rj)} with the func- 
tional derivative to be interpreted as a directional derivative along the surface of paths 7. The 
longitudinal gauge variables have a complicated formal implicit expression given in Eq.(4- 
49) of the second paper of Ref. |TJ and satisfy {t]^'(t, a), T b (r, a)} = —5 ab 5 3 {a — a), where 
r o (r, a) = T b (r, a)A ba (r]^ A \T, a)) are the Abelianized Gauss laws [{T a (r, a),T b (r, a')} = 0]. 
In Eq.([73|), A ab (j]^ A \T,a))dr]l (r,a) is the pure gauge part (saturated with da it is the 
BRST ghost) of the vector potential A a (r, a): the magnetic field B a (r, a) is generated only 
by the second term of Eq.(^). In this sense, % '(t, a) = is the true generalized non- 
Abelian Coulomb gauge with all the same properties of the Abelian Coulomb gauge. In 
suitable weighted Sobolev spaces, as discussed in Ref. 0, this gauge-fixing is well defined, 
since all the connections over the principal SU(3)-bundle are completely irreducible [their 
holonomy bundles (i.e. the set of points of P(R 3 , SU(3)) which can be joined by horizontal 
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curves) coincide with the principal bundle itself] and there is no form of Gribov ambiguity 
(i.e. of stability subgroups of the group of gauge transformations for special connections 
and/or field strengths). In these spaces, the covariant divergence is an elliptic operator 
without zero modes ]nj and its Green function Cab\3, 3'; r) is globally defined 

A (A) 

D ab (r, 3) ■ &\a, a';r) = -5 ac 5 3 (3 - a') 



a'; r) = c{3 - a'; r) = c(3 - <?){PeP # ) ab , 

c{x)=dc{x) = ^{x) = ^, A = -2», c(x) = y%x) = ^L. (74) 

The path ordering is along the straigthline (flat geodesic) joining y and x when S T is a 
hyperplane like £w>- 
Therefore, we have 

d - 

vr a (r, a) = -—d ■ 7? a (r, 3) + n a± (r, a) = 

r N 
= n a , D± (r,3) + I d 3 a ^i\a,a'-T)[T b (T,a)-J2Qib(r)S\3' - ff^r))) 

i=l 

d ■ TT a± (T, 3) = D ab (r, 3) ■ n biD± (r, 3) = 0. (75) 
It is shown in Eqs. (5-7), (5-8), (5-10), of the second paper in Ref. that we have 
d-7r a (T,a) = J d 3 a C^\a,a;r) 

N 

[c bef A e (T, a) ■ v?/±(r, a) + T b (r, a) -^2Q ib (r)5 3 (3' - ^(r))], 

i=i 

3) = J d 3 a [S^5 ab 5 3 {3-3')- 



" 2?(*, r)c ceb Ai(r, 3'))4Ar, <?') 



^KAr^)=P l l^K D ±(r,3), (76) 

where P±(3) = 8^ — <9 4 <9 J /A. Moreover, Eqs. (5-21) and (5-25) of that paper give 
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R,D±(r,a),r 6 (r,a')} = 0. (77) 

Therefore, the color SU(3) canonical pairs of Dirac's observables turn out to be A a ±(r, a), 
Tt a ±.(j ■,<?)■ They satisfy the Poisson brackets 

{*Ll(t, a), fl ± (r, a')} = -6*1* ~ <?). (78) 

Instead, the gauge sector is given by the pairs A aT (r, a), 7r^(r, a) ~ 0, t^{t, a), r a (r, a) ~ 0. 

Let us now look for the Dirac observables of the particles. First of all, the Grassmann 
variables are not gauge invariant because one has 

{9 ia (r),T a (ra)} = (T a ) a(3 9^5 3 (a - ff^r)), 

{d* a (r), Ura)} = -0*^)^6* (a - #(r)). (79) 
The Grassmann Dirac observables are 

^(r) = [Pe^ ) ^ ) ^«)^]^(r), 

=>{e ia {T),r a (T,a)} = o, 

&(r)=fl5,(r)[Pe n 2 ) (^ ) (-.-«W))^] /Ja , 

{e ia (T),e* jP (T)} = -iS ij s aP , (so) 

with the path ordering evaluated in the fundamental representation of SU(3). The Dirac 
observables for the non-Abelian charges of the particles are 
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Qia — i9*a{T a ) a i39ii3 — 

= Q ib [Pe n ^ ){v{A>{T ^ {T)))fa ] 



{Qia,r b (r, a)} = 0, 

{Qiai Qjb} $ijCabcQia (®-Q 

where in the second line the path ordering is evaluated in the adjoint representation of 
SU(3), since one has used the identity 

e u b T» T a e u cT c = T c ( e -^) ca . ( 82 ) 

By using Eqs.(|75]) and ($[§, for ^ A) (r,a) = T a (r,a) = [so that also dr)i A \r,a) = 0, 
(?/ A ) (r, ex) ) = as shown in the second paper of Ref. namely in the generalized 
Coulomb gauge [we still go on to use the notation {.,.} for the new Dirac brackets with 
respect to the second class constraints rj^(r,a) ~ 0, T a (r,a) ~ 0; instead the temporal 
variables A ar and n T a simply decouple], we get 

vf a (r, a) -> 7r a (r, a) = n a± (r, a) - 

a N 

^ J i=i 

Cl^(<?, *;r) = c(a - a')(Pel> *^*&)*. (83) 

While in the electromagnetic case it is possible to get the physical Hamiltonian with- 
out imposing the Coulomb gauge- fixing fj em (x) ~ as in Ref. |F| (namely it is obtained 
by a canonical decoupling of the gauge degrees of freedom ), this is too difficult in the 
non-Abelian case. Therefore, we shall evaluate the physical quantities by imposing the gen- 
eralized Coulomb gauge-fixings rjf^' (r, a) ~ 0. Conceptually, the canonical decoupling of 
the gauge degrees of freedom gives the same results for the physical quantities. 

As shown in Refs. ]]],§], the Noether identities implied by the second Noether theorem, 
applied to the color SU(3) gauge group, give the following result for the weak improper 
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conserved non-Abelian Noether charges Q a and for the strong improper conserved ones 

Q(s)a 

N 

Qa = g- 2 c abc / d 3 aF b ok (r,a)A k c (r,a) - Y,Qia= 
J i=i 

= Q(s )a = J d 2 Z-E a (r,a), (84) 

Then, we get [see Eqs.(6-33)-(6-35) in the second paper of Ref. M\ the following Dirac's 
observables 

N 



i=l 

N 

= / d 3 ap a (r,a) = / d 3 a[Pa YM) (r,a)+Y,PU^^ 

J J i=l 

pi YM) ( r > = -c afec A 6 _L(r, a) • 7f c _L(r, a), 
Pia( r > <?) = ~Qia(r)5 3 (a - ffi(r)), 

{Q a , Q b } = c abc Q c , {Q? m) (t), Q { b YM \r)} = c abc Q? M \r), 
{A aJ _(r, a), Q b } = c abc A cJ _(T, a) 

{v?a±(T, a), Q b } = c abc ^c±(r, a). (85) 

While particle's 3-positions ffi(r) are gauge invariant, for particle's momenta we have 

K(r), T a (r, a)} = -Q m ( T )^6 3 (a - ff t (r)). (86) 

To find the Dirac observables Ki(r) [namely to dress the particles with a color cloud], 
we note that, anaguously to the electromagnetic case 0, we expect that the form of the 
minimal coupling is mantained, namely 

«i(r) + Q ia (r)A a (r, = ^(r) + Q ia (r)A a± (r, ^(r)). (87) 

Since we have 
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Qi^A^Ur)) = Qa^^e^^'^^^^l^il^^r)), (88) 
Eqs.([73|) suggest the following expression for the gauge invariant momenta 

ki(r) = Kiir) +^g, a (r)e a (7 7 ^(T,77 i (T)),^(T,^(r))). (89) 

a 

Since we have 

{e:(r ] ^\r,a),dr ] ^(r,a)),T b (r,a')} = = -B uh ^ A \r, a')) ^f^^ *) = 

f)r]« ; (r, a') 

= c abc e r c (rj^(r, a), drf A \r, a))5 3 (a - a') + 

+ S <* Q^r ' ( 9 °) 

we get the gauge invariance of /^(r) 



K(r),r b (r, ( x)} = K(r),r b (r ) a)} + 

+ {Q ta (r)Q^ A \rMr))M A \rMr)))Mr^)} = 
= -Q l b(r)^5 3 (a~ff i (r)) + 

+ c abc e r a (^ A \r, fu(T)), d V ( A \r, ^(rmUT^ia - iftr)) + 

+Q J a(r)c abc e^(r ? ( A )(r,r/,(r)),^ A )(r,r7Kr)))5 3 (r/ J (r)-a) + 
+ Q lb (r)^5(ff l (r)-a) = 0. (91) 

We can now rewrite the constraints 7"^(t) « 0, 7^ p (r) « 0, « in terms of the Dirac 

observables with respect to the SU(3) gauge transformations. The Grassmann constraints 
become 

^ = Efc(*(r)«0. (92) 

a=l 

As shown in Ref. QXJ] , we have J2 a Bl(r, a) = J2 a B a (r, a) with the chromomagnetic field 
B a (r, a) built in terms of A a j_(T, a). For the chromoelectric field we have from Eqs. (|83|) and 
using Eqs.flgoD and ([F|) 
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E <{t,*) 



r,[ A) =r b =o 



= a) + | J So A • 2^>(*, * l5 r)p b (r, * x ) f = 

= E[^-l( t > s ) - ^ / d 3 <7i(5 ab <5 3 (<7 - <Ti) + c amn A m± (T, a) • Clfe x) (^ ? r ))Pf>( r > <?i) ] 2 > 

"'a 6 •'a 

+ / d 3 aYll^J d3 °l( S «^-Zl)+CarnnArn±(T,B)-^ (93) 

To obtain the last line, we have done an integration by parts and used the transversality 
of 7r a ±(T, The final result is 



4 A) =r 6 =o 



S s 2 /^ E^(r^) 

J a 

= f d 3 aJ2€±(r^) + V[ff i ,A a± ,fc a± }(r), (94) 

J a 

with the potential given by 

V[rji, Xa±, fa±](r) = g 2 s J d 3 ad 3 a,d 3 a 2 ]T [d a ^\ff, *x\ r)A,(r, a,)] 



9 2 s j d 3 a 1 d 3 a 2 E Pa(r, <?i) K^fii, di\ t) p b (r, a 2 ) = 

a,b 

N 

g] / fa^YX^h) "E&aOWi - #(r))] 

J a,6 i=l 

V 

^ab(<?i, ^2; r) [pi yM) ( r > 30 - E Oj&W 3 ^ - %(r))], 



j'=i 



Kab(a 1: a 2 ; r) = X 6a (<7 2 , <?i; r) = 

<U<5 3 (<7 3 - <?i)5 3 (<7 4 - 02) 



= y rf 3 <j 3 (i 3 (T4{ 



47T I (7 3 — O4 



+ 



^(^-^^(r.^-C^fe^jr)]^ 
H : — —\ h (0-1 ^ (T 2 J + 

47T I (7 3 — CT4 I 

A.(r, ff 3 ) ■ C^fo, gij r)Ufe gj) ■ g^jgj, <? 2 ; r)k , Q _. 

+ -r-^ — • ( 95 > 

An I cr 3 — CT4 j 
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The first line of this equation agrees with Eq.(6-25) of the second paper in Refs. |T|, while, 
from Eq.(6-27) of that paper, we get that K ab (3i, 3 2 ;t) = —G^^(ai, 3 2 ;t), where G^ A ^ is 
the Green function of D ac (r, 3) ■ D cb (r, 3) [see Eqs.(3-16), (3-20), (3-21), (3-25) of that 
paper] . 

The constraint giving the invariant mass of the system takes the form 



U{t) = e s - Vi\™i + fair) + ^ Q^{r)A a± (r, 7fc(r))) 



2 



^2 



rJ2 



(X (X s 

= e s - H rel « 0. (96) 

We see that, since p a = p^ M ^ + X^i Pia, there is a universal interaction kernel 
Kab(&i, ®2\ t) which creates the particle-particle, the particle- field and the field-field in- 
teraction between the corresponding color charge densities. This interaction kernel contains 
3 kinds of instantaneous (i.e. at equal r on the Wigner hyperplane) interactions: 

i) a Coulomb interaction; 

ii) an interaction mediated by an arbitrary center (over whose spatial location is in- 
tegrated) : one color density has a Coulomb interaction with the transverse potential at 
the center, which simultaneously interacts with the other color density through an instan- 
taneous "Wilson line" along the geodesic (the straightline) on the Wigner hyperplane, i.e. 
[Ii = ~A a i_f a with (f % = c cab }: 

[1l(t, a) ■ (^(3, 3; r)] ab = A c± (r, a)c cad ■ c0 - 3) (p e />'^^'^ ) db , 
where a is the position of the center and a that of the color density; 

iii) an interaction mediated by two arbitrary centers (over whose spatial location is 
integrated): each color density interacts, through a Wilson line along the geodesic, with the 
transverse potential at one center and the two centers have a mutual Coulomb interaction. 

If we rescale the transverse potentials and electric fields [A a ± = g s A a ±, 7r a ± = g^TTaX, 

=^ P { I M) = ~ C abcAb± ■ 7T C _L = C abc A b± ■ 7T C± ], We get 



\ 



mf + {Ki{r) +^^Q a i(r)A a± (r,^(r))) 5 
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E / £ *) + ^ " 9s ^ 97^a±}(r) « 0. (97) 



a a 



A 



To get the constraint defining the intrinsic rest frame, i.e. 1 H p {t) ~ of Eq.(^), we must 
consider the term J2 a ^aij, B) x -B a (r, a) = J2 a ^aiji B) x 5 a (r, <?). From Eq. (|83D we get 

y d 3 <r£ 7r a (r,a) x B a (r, a)|^) =rfc=0 = 

= / ^^[^ a± (r,a) + - j • ^(^aijr)^^,^)] X 5«(r,5). (98) 

Remembering that (7? a x B a ) r . = e rmn ir™B% = F ars Tr s a , we get 

J d 3 a K{r, a) x B a (r, B)] r \ v (A) =Vb=0 = 

= f d 3 aJ2Fars(T,a)n s a± (T,a) + 

J a 

+ / d 3 <7^[<9 r A a _u ~ d s A a±r + c abc A b±r A c±s }(T, a) ■ 

J d 3 aid a ■ Cld x) (^ ar, r)p d (r, B x ) = 
= [ d 3 aJ2F a ±rs(r,a)r a± (r,a) + 

J a 

+ / d 3 aJ2Ab±r(T J B)[5 ba A - c bca A c±s (T, B)d s ] ■ 
■ 1 J d 3 aA ■ ax; r)p d (r, a,). (99) 

We have done a first integration by parts to use the transversality of A a ± and then a 
second one on d s A a j_ r . In the last expression we recognize the Faddeev-Popov operator 
■K&a (r, &) = -d- D ha (r, a) = 5 6a A + c bac A c± (r, B) ■ d = 5 ba A - c bca A c± (r, 3) ■ d. As 
shown in Eqs.(3-ll) and (3-17) of the second paper in Refs. 0, we have k ( £ ± \T,B)-^ r d a ■ 
Qc(B, B\\ r) = -5 ac 8 3 (B -B x ). Therefore, we get 

J d 3 a J2 [*a(r, B) x B a (r, B)\ r |^A) =r&=0 = 

= / d 3 aJ2[FarsK± - A a ± r pa] (t, B) = 
J a 

= J ^^[(d r A a ± s — d s A a ± r + c abc A b ± r A c ± S/ 



7T 



a± 
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N 

A a ±r(-Cab c A b ± ■ 7T C _L + ^Pia)]{r, B) = 

i=l 

N r 
= EEQ-(^±r(T,^(r)) + / d 3 aJ2[(d r A a±s -d s A a±r )r a± )(r,a). (100) 

i=l a J a 

Then, the constraints H p (r) ~ have the form 

N r 

W P r(T) = EM T ) + / d 3 aY,[(d r A a±s - d s A a±r )r a± }(T, a) = 

= k+rif) + I d 3 aY,[(d r A a±s - d s A a±r )7T s a± ](r,a) « 0. (101) 

As in the electromagnetic case of Ref. 0, there is no interaction term in this constraint 
(neither minimal coupling to particles nor self-coupling of the color field). As shown there, 
this is the requirement for being in an instant form of the dynamics. The other requirement 
is that the rest-frame spin tensor 5J S must also not depend on the interaction. From Eq. ([71]) , 
rewritten in terms of ff^ Ki, one has 

S rs (r) = J2 VM + E^E QiaAim - (r <- a) + 

i i a 

+ Jd 3 aa r F r au (a)7r:(a)-(r^s) = 

= E Vi*l + E^E &Al (ffi) - (r <- a) + 

i i a 

+ J d 3 o a r F s au {a)*l{a) - (r <- s). (102) 

By doing various integrations by part, discarding terms symmetric in (r <->• s) and using 
the Faddeev-Popov operator, we obtain 

J d 3 a a r F asu (a)fc^(a) - (r s) = 
= J d 3 a a r F asu (a)^ ± (a) + 

+ J d 3 a a r F asu (a)^- J d 3 a' d ■ ( ab (a, a')p b (a') - (r s) = 
= J d 3 a a r F asu (a)jr^ ± (a) + 

- J d 3 a d u (F asu a r ) J d 3 a' 9 ' {<?) - (r ~ s) = 

= J d 3 a a r F asu (a)n"Z ± (a) + 
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— I d 3 a [-d,A arl (3) - d r A all (S) - d u d„(A all (a)a') + c ahk A k ^(3)d'(A hll (Sy)} ■ 
= J d 3 a a r F asu {a)n u a± {a) + 

+ J d 3 a J d 3 a' [d u d u (A as± (a)a r ) - c hka A ku± (a)d u (A hs± (a)a r )} 9 ' Cab ^ 7 ^ p b (a') + 

— (r <-> s) — 

= J d 3 a a r F asu (a)^ ± (a) + 

— # — # 

+ y><7 |dV A hs± a r (a)[d u d u S ha + c hka A ku± ([cr)d u ] - ' Cab ^ ^ p b {B') + 

— (r <-> s) = 

= y d 3 <r ^^(^^(ct) - J d 3 a A hs± (a)a r p h (a) - (r <-> s). (103) 
In conclusion, we get [in the last line we use H p (t) ~ and go to relative variables] 
^' = E^K( 1 -)+ / d 3 aa r J2id s A a±u (r,a)-d u A s a± (r^)K±(r,a)- 



= 1 

(r <-> s) ~ 

AT-l 



E A«(r)ifS(r) + / d 8 ^ - i£(r)) Et^A^r, a) - (U^r, <7)]t&_(t, a) - 

a=l J a 

(r <-> s), (104) 



expected in an instant form also in the non-Abelian case. 
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V. THE REDUCED HAMILTON-DIRAC EQUATIONS 



To write the reduced Hamilton equations, it is convenient to introduce the gauge-fixing 

X = T s -r^0, (105) 



whose conservation in time requires A(r) = —1 in Eq . (|72|) . Then, we can eliminate the 
pair of variables (T s ,e s ) and describe the evolution in terms of the rest-frame time. The 
Hamiltonian for this evolution will be 

N 

H D = H rel - A(r) • H p {r) + ]T ^(r)iV 4 , (106) 

i=l 

with 



Hrel = ZW 771 * + («<( r ) + Z Qia{r)A al (r, ^(t))) 2 + 
i V a 

+ Z^J d a \ 2 2^2 J + 2 ^' ' ^ ^ ^ 

If we would know the correct form of the gauge-fixing to eliminate the intrinsic center-of- 
mass degrees of freedom associated with the constraints 'H p {j) ps 0, we could put A(r) = 
and rewrite the invariant mass H re i only in terms of relative variables. See the analogous 
discussion for the electromagnetic case in Ref. |L4| 



At this stage of reduction, we get the following Hamilton equations of motion [r = T s ; 
PJ_ s (a) = 5 rs + d r d s / A; ffi = (rjl); d/dffi = (d/drjl); the symbol = means evaluated on the 
solutions of the equations of motion] 

-rVi{r) = r)i , , = = - A(r) 

dT V™? + fair) + EaQ ia (r)A a± (T,ff l (rW 

^«*W = -E(-^ + ^ W)EQ*W ^ 

-%2^W r ) / "q^ Pb{T,a'), 

a,b J OT li 
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- g 2 s I d 3 aK ab (ffi(T), a; r)p 6 (r, a)] - i/ii(r)9i a (r), 

3* (t) 

dr 



g s / d aK ab (ffi(T),a;T)p b (T,a)]+iHi(T)9* a (T) 



-j-Qia(r) = c acd Q id (r)[^(r) ■ A c ±(T,ffi(r)) - 

CLT 

-9s J d 3 aK cb (fj l (T),a; r)p b (r, a)], 

(108) 



— A a _L r (r, a) = -g s TT a ± r (r, a) - [A(r) • <9]y4 a _L r (r, a) + 

+ ^ 2 Pl s (a) | ^ff'catatiij.Cr, 5)^(5, ^;r))]p c (r,a') 
|l^±(r, *) = ^ 2 Pr(a)^ afc (i ± )(r, ^(r, a) - [A(r) • ^(r, a) - 
- E ^a(r)Pr(a)(^i + A s (r))5 3 (a - #( T )) + 

1 nfs ^ /^K» A±,7Ta±](r) 

I oA a±s (r, a) 

d d N N 

f Qi™)(r) = = -c acd Y,QUrm.(T)-A c± (TMr)) ~ 

-2 2 | rfV^ cb (r/,(r),a;r)p fe (r,a)], (109) 

since J^Q a = [it is not known (see also Ref. [[[]]) how to check this formula by a direct 
calculation]. 

In Eqs.( [109| ) we have 

lSV[fJi,A a± ,7l a± }(T) 2r>rs/-\ ~s I —*\ fj3 ' rs t- -> \~/ ->'\ , 

o Tl 7 V = -g s PKv)cabv< ± {T, a) / d 3 a K bc {a, a ; r)p c (r, a ) + 

2 M a _L S (r,(r) - - J 

+9 2 s PT(v) J d?a l d 3 a 2 p b (T, ai) 

{ / , ,f a ., [-^ - ^)c a6e Ci^ )s (^, ft; r) + 

J 4TC\a — (72 1 

+c be/ 2 /i (r, a) ■ , ft r)c man CH ±)s (a, ft; r)] + (ft <-> a 2 ) + 

+ / J,^ 4 ■ - W^'fo, ft; r) + 
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+c bef A f± (r, ff 3 ) • Qi x \a 3 , r)c man C^ ±)s (c?, a i; r)] • 
c e r S ^ s ±(r, a 4 ) ■ C™ x) ((?4, ^2; r) + 
c bef A f± (T, a 3 ) ■ ^ x) (a 3 , #\\ r) 
\-5 3 (a - a 4 )c acm Cil x)s (^4, o 2 - t) + 

+C cmn A n± (T, a 4 ) • Cmr X) (^ °\ T ) c ras(ii x)S (v, °2] t)]} }/5 c (t, ff 2 ), 

(110) 

because from Eg. ([Tip we get 

dA u±s (T,a) 

^ —Jl < ~\ = P ?( a Kbe {<7i,<7;Vc eud Q c x> (a,a 2 ;r). (Ill) 

These equations of motion have to be supplemented with the constraints 



N 

n nr (r) = Y, Kir(r) + / d 3 aY,[{d r A a ± s - d s A a±r )r aA ]{r, a) « 0, 



L p r V 

i=l 



a=l 



The first line of Eq. (|108 ) can been inverted to get 



fti(r) ^tpm . ^ (r)+A(r J -J2QUr)Iax(r,m(r))- (113) 
v/i - (m(T) + Mr)) 2 - 

The first line of Eqs.( |109| ) coincides with Eq.(6-24) of the second paper in Refs. |J 

if in this equation if a ± — > — 7r a j_, iip^T c ip — > Qic an d if we note (see Ref. that 

D^ x)k {r t a)P^{^ x)n {^)r) = 0, K^uS^t) = -G^{ff 1} ff 2] r) and (see Eq.(3- 

17) of that paper) 3^>(<?, <?'; r) = -^(r.^G^^^r) = i£/V. Wc(*,*V). 

Therefore, Eq.(6-21) of that paper gives its inversion in the form 

K±(r,3)=g; 2 E: ± (r,S)= - 9?P?(3) f d^V^'US,^; r) • 
[^ + A».A]i. ± (x, 3l) - 
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N 



-Pl s {a) J dV^^^r)^^) = 

8=1 

= -g: 2 Pl s {a) J d 3 ax[5 st S ab S 3 (a-ax) + 

+ D { a i- )s (r, a)K de (a, ft; t)D^\t, ax)] + A(r) • + 

+ Pri^CadcA^^KMM^r^ir), (114) 

where [see Eq.(3-30) of 

Ref. i] V { ^ )rs {B^r) = 5™5 ab 5 3 (a - ax) - ^ )r (r,a)G^(a, a i; r)^ 5 ^^) = 
5 rs 5abfi 3 (& — o\) + Cac" L ^ r (^ T )Dit ±S)S ( T > ^i) * s a projector giving an explicit realization 
of the Mitter-Viallet abstract metric [it satisfies / d 3 axV ab ^(a, ax, T)V bc """^(ctl, tr'; r) = 
PfrWfr^T); Id^xV^^a^x^D^^ax) = I ^(r,^ (<?i,<?;r) = 
0]. 

By using also Eq.(6-23) of Ref. [|, namely £>if ±)r (r, a)P r ^{a)b { ^ 1 - )s {t, a) = (when 
acting on functions of a), we get 

p a (r, a) = & Wfcifr ^(a) | dV^^a, ft; r) [^ + A(r) • Jr]iL.fo ft) + 
iV . 

+ E / dV^C* - ft) + c abc A r b± (r, ff)Pr(Z)( { J x)s (<?, ft, r))p id (r, ax) = 

i=l 

_ d - d - N 

= g s 2 c abc A b± (r, a)-[— + A(r) • — }A c± (t, a) + ^/5 ia (r,ft) = 

N 
i=l 



p a YM \r, a) = -g; 2 c abc A b± (r, a) ■ E c± (r, a). (115) 



Therefore by putting A(r) = g(r) like in Ref. |L4| , the equations of motion for A a ±(r, a) 
become 

d , d 



+ Mr) ■ ^}{-g7 2 Pl s (f) J d 3 axVH ±)st (t, ft; r) • 

[-^ + t(r)~}At ± (T,ax)- 

r - N 

-P?(a) J d 3 ax(^\a,ax-,T)j:p ib (T,ax)}~ 

i=i 
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-g- 2 P^{a)[5 ab A + c abc A c ±(r, a) ■ d]A s b± (r, a) == 

N 

= E Q ia (r)Prmvt(r) + \ s (r)]5\a - ff^r)) + 
i=i 

-g 2 s PY{^)c a bvK±{ T ^ v) j d 3 a'K bc (a, a;r)p c (r, a) + 
+g 2 s P]_ s {a) J d^d'a^r,^) 

{ I J ,t a ., [-^ " *)c abe (^ )s (v, Bi\ r) + 

J 47T|(T — 02 1 

+c be/ 2 /± (r, a') • C^f/, *5 r)c man C^ )s (^, <? l5 r)] + (a, <-> a 2 ) + 

+ J 4 t\?f-% {[ - 63{ff ~ ^Wi^'fo, ft; r) + 

+C bef A f± (T, <7 3 ) ■ Cim^^S, ft r )CmanC^ ±)S (0 ; ) ft; 7")] • 

c ers A s± (T, ft) ■ C ( ^ x) (^4, ft; r) + 
c bef A f± {r, ft) ■ <f^ x) (ft> ft; r) 
[-5 3 (a - ft)c acm Cm U ±)s (ft, ft; r) + 

) • Ci A /H^4, ft r)c ras Ci^ s (ft a 2 ; r)]} }p c (r, ft), (116) 

After some manipulations its final form is 

P ?(v) + t(r) ■ ^] / ^^'(ft r)[^ + g(r) ■ A] + 

= - <#T(^{E&.(r)[ij?(r) +r(r)]5 3 (a-^(r)) + 

i=l 

+ c afed irf_L(r, a) /" d 3 a'K bc (a, a ; r)p c (r, a')} - 

- 9iP±{*)AdL(T, a) J d 3 a 1 d 3 a 2 p b {r, ft[F abc(i (ft ft, ft; r) + 

+ CauvC v d e K eb (cx,ai;T)K uc (a, ft; t)]/5 c (t, ft), 



or 



Pl s (t) {[§p + iT(r) • ^] / d^^V - 5) + 
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+ D^{r,a)K uv {a,B-r)b^\r^][^ + fa) ■ ^] + 
+ 5 st J d 3 a5 3 (a - B) [5 ad A + c ad J c± (r, B) • ^] }A t d± {r, B) ± 

N 



i=\ 

+ CaM^dj.(r, a) J d 3 aK bc (a, a ; r)p c (r, a)} - 

- gtPl s {v)A s dl _{r, a) j d 3 a 1 d 3 a 2 p b (r, ai)[F afccd (a, a 1 ,a 2 ;r) + 

+ c auv c vde K eb (a ■,&i;t)K uc {o ', a 2 ; r)]p c (r, a 2 ), 



or 



P r ± s (a) {S st [6 ad ((^ + fa) • + A) + c ad J c ±(r, 

J d 3 a5 3 (a + + fa) ■ ^\D a tfa, a) J d 3 a 
K uv (a, B; t)D&»(t, B)}[^- + fa) ■ JL] }A a± (r, B) ± 

N 

° „2 r>rs 



~ 9 2 sP1 s (°){EQUtM(t) + 9 s (r)]5 3 (a - ff^r)) + 
+c a bdA d± (r, a) J d 3 a'K bc (a, a;r) 

N 

[pf M) (r,a') + EA C (^-)]}- 
i=i 

r N 
-9lPl S (v)A d± (r,v) / d 3 a l d 3 a 2 [p b YM \r,a l ) + Y,p lb {r,a l )} 
J i=i 

[F a bcd((?,(Ti, <?2; t) + c auv c vde K eb (a,ai;T)K uc (a, a 2 ;r)] 

N 

[pf M) (r,a 2 ) + Ep, c (^^)], 

i=i 



TP /- - - \ ( K fc {a,&i,T) K fc (a,a 2 ;r) 

rabcd\V,Vl, °2\T) = -C abe C ed f{ — - + —-) + 

4:71 \(T — a 2 \ 47T|<7 — <7i| 

dV 

+ ' 



( K kc (a, ai;r) K kc (a,a 2 ;r) 
47r|(7 — cr 2 | 47r|cr — <7i| 
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- / — -rz — —M abe c edf K fu (a, fix] r)c ers A s _L(r, a) ■ Q w -^(<t , <? 2 ; r) + 

J 47T|(T — er | ' 

+ c be /A/_L(r, a) ■ Q£ x \&, &i] r)c acm c mdn K nu (a, a 2 ] r)] + 
/* tZ^ ^"2 v j j 

+ y 4?r | (? 1 _ J 1 [c6e/^/±(T, ^3) • Qto X) (^3, T)c ma nCndkKku(v, T) 
C ers A s± (T, B4) • Cru X \^ °l\ t) + 

+ c bef A f± (T, a 3 ) ■ Ci^ x) (^3, r) 

C cm „i«l(r, CT4) • <?; r)c ras C sdk K ku (a, Ol\ t)\, (117) 

where the identity P r ± s (a)( a i ±)s (a, a'; r) Pl s (a)£>(^) s (r, 5)^(5, a; r) 

= PX s (a)c acfe ^_ L (r(T)i^ C 6((j, a; r) was used in Eq. (|110|) to obtain F abcd . 
The equations of motion for the particles are 



dr 



d farm i -^Q fa (r)4(r,jji(r))]^ 
= " mr) + <f(r)] E g to (r) ^^ (T)) - 

a a 7« 

9s2^Qia(r) J d a a — Pb(r,cr), 



a.b 



or 



A farm . !!M±M4 } ± ^ l a± ( r , 77,(r)) Cacd 4 d (r) 



K( r ) ■ rfi(r)) - g s J d aK^^r), a; r)p b (r, a)) + 
+ £<?*(r)[ Wr + %(r) m ] - 

2v-n I \ f j3 dK ab (ffi(T),a; t) . 
-# s E^( r ) _/ d a ^ Pb[T,a), 



or- 



= EQM{KAr,Vi(r)) + mr) + 'g{r)) X B a {r, for))]*} - 

a 

- E Q l a(T)Pl s (fJ i )c amd A s m Ar, f)i(r)) 

a 

J d 3 aK de (ffi(T), a; T)c enb A n± (T, a) ■ E b± (r, a) + 

N 

a,b j=l 

+ aJf - tW( I ) ^ (r)ir) fe(r) - 

-Qia(r) J ^alCadcA^^ffiiT^K^fJiiT),^; t) + 

dK ab (ffi(T),a;T) {YM ) { ^ 
+ ^ }Pb (r,a)}- 

-J2QUr)'g(r) ■ [ gjkfe^jj _ Cac6 l c± ( T ,^(T))i^(T,^(r))] }, (118) 

a °Vi 

where Eqs. (|108D , (|T09D , (TO) , QITSD , the definition F a rs = 9 r i^ ± - <9 s i^ ± + c a6c I^i^ ± = 
e rsi S* and D a j ±)s (r, a)Pf(a)l)^ ±)r (r, a) = have been used. Let us remark that Eqs.flIT7D 



and (|118|) are a system of integrodifferential equations due to the presence of the transverse 



projectors and a completely open problem is how to define an initial data problem for them. 

Eq.( |118| ) is the non-Abelian version of the particle equations of motion given in Ref. 
fl4| , which should produce the Abraham-Lorentz-Dirac equation if we were able to get a 
non-Abelian Lienard-Wiechert potential from Eq. (|117|) . The last term would be absent in 
the gauge A(r) = g{r) = 0. This equation shows explicitely the chromo-electric [E a ±] and 
chromo- magnetic [B a , A a ±] forces acting on the particles. Moreover, the explicit interparticle 
instantaneous color forces are given by 

4a(T)E^(r)[c^Jr,^(r))ir c6 (^(r),%(r);r) + g^^(lM ] = 

a,b j=l 0r li 

N 

a,b j=l 



[cadcA d± (T,fji{T)) / dV 3 d J <r 4 { 



Att I <T 3 — a 4 
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+ 4n\S 3 -fft\ +Vk»r,j) + 



<9??[ J " ' 47T I <7 3 — <T 4 

£ 3 (a 4 - ^(t))[1 ± (t, ff 3 ) • ft^fo, fH(r); r)] cb 
+ 4vr | <r 3 - <? 4 | + {Vl Vj) + 

47T | (T3 — (74 | 

It is evident that there are divergences due to the self-energies for i=j. They will be 
discussed in the case N=2 in the next Section, because their regularization is connected 
with the existence of asymptotic freedom already at this pseudoclassical level. 

After an integration by parts, the constraints H pr (T) « become 

Efam, . ^ T)+ ^ T) - £Q^A±(r,^(r))] + 
\A - (^(r) + g(r)) 2 * 

+ g; 2 f £a[(dA s a± )E s a± ](r,a) = 0, (120) 
and the expression of the constant invariant mass H re i becomes 

N 

x . riirrii 



2 



i=l \/l-(f}i(T)+t(r)) 
+ I^,X^la±]+E^I / ^MC+^IK^)- (121) 
To find the classical theory implied by the pseudoclassical one we follow Ref. |15[], where 



the Berezin-Marinov distribution function [see in Refs. [0] for the Grassmann variables of 



phase space is defined for the case of SU(3) [in this discussion we will denote 9 the Dirac 
observables 9]. For each particle i, on the space of analytic functions of the #; Q 's, namely of 
the functions f{6i) = a% + PiaOia + lia^-ap^OipO^ + Sie a /3 7 9i a 9ii39i 7 , the more general density 
function pi(9i,9*) satisfying the normalization condition / pidpi = 1 [dpi = Ila=i d9 ia d9* a is 
the Berezin integration measure on the Grassmann sector of phase space] and the positivity 
condition / p,J* fdp. t > for any analytic function f, is [q = q(t), Wi a p = Wi a p{r), Vi a p = 
V ia p{r)] 
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Pl (6 t , 9*) =c l + Y: e: a W iaP 9 if} + 2(Tr Vi) (£ 9*J ia ) 2 - 4(£ 0* a V ia(3 9 ip ) (£ 9*^) + 

a/3 a a/3 7 

+ fi(XX*a<^) 3 , 
a 

Ci = c*> 0, 

Wi = W} = w io lL + w ia -X a , definite positive, 
Vi = V} = Vi H + Vi a -X a , definite positive, 



w io =^TrWi, w ia = Tr{\ a W l ), 

v lo = ^TrV l , v ia = Tr{\ a V), (122) 

since one gets / pif* fdpi = \a\ 2 + 36ci\5\ 2 + A(3*Vi(3 + A^Wi^. Other useful formulas are 
[\d a bc = Tr [{|A a , |A b }|A c ] ] 

/ Pi^iJifsdpi = AVif3 a , 

J Pi9* a 9i l3 9* J 9isdpi = Wi\ T e ai \6ps T , 
J Pi9i a 9if 3 9* J 9is9* x 9i T dpi = Cie ai \6ps T , 



< JVi > = < XX0i« >= 12 ^o, 

OL 

< Qia > = < X]^*a(o^o)"/3^ >= ^ia, 

a/3 1 

<N 2 >=< (E««) 2 >=Qw lo , 

a 

< (£C(;U a )a/3M(£ W >= -J^a, 

a/3 Z 7 Z 

<iVf > = < (E«a) 3 >=6Ci, 



< QiaQib =< ((^6* a (^K) a p9ip)(^9* 1 (-\ b ) 1 s9is) > 

a/3 ~" -yS 

= \{.d abc w ic - w io 5 ab ), 
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< QiaQibQic >-< (^Oi a (^K)ap8ip)(^0i^h)jS0is)(^6i^K)fiu8i V ) >~ 

< Q ia Q lb N t >=< E^)^E^4^WEW >= 

= -^(W (123) 

Therefore, the constraints iVj w imply at the classical level [at the quantum level in 
general q^O due to ordering problems] 

Ci = v io = w io = w ia = 0, 



PiiOiJD = -AViaWQiatfNi + ^N? « 0, 



^ Qia ^* 2fj a , < QiaQib ^* ^ QiaQibQic ^* 0, 

< JVi >=< JV? >=< JV* >=< QiaQ^JVi >= 0. (124) 
The classical observable associated to an even Grassmann-valued function g is 

- N N N 

<9>= g~[[pidpi = 9pdp, P=Y[pi, dfi = Y[dfii. (125) 

J i=l J i=l i=l 

The distribution function p = n£Li pi must satisfy the Liouville equation 
^ + {p,H D } = 0, 



Pl-Pfc-l[^ + {Pfe> ^£»}]Pfc+l-PiV = o, 

fc=l ° T 
{p , M p (r)} » 0, 

{p , iVi = * 0, i = l,.,iV, (126) 
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with Hd the Hamiltonian of Eq. (|106|) . The last two equations are identically satisfied. 
Therefore, by using Eqs.( |122j ) and ( |108| ) for Q ia , we get the equations [valid in a neighbour- 
hood of the region iVj « 0, where however we will put Nf = in the sense of Dirac's strong 
equality] 

N 

Y,Pi-Pk-i[-4{V ka f>+ [V k ,T a } a/3 [ff k (r) ■ A a± {T,ff k {r)) - 

k=l 



-9s J d 3 aK ab (rf k (T),a; r)p[ yM) (r, a) + 



N 



+g 2 s £ K ah {rj k {T)Mr)\ r)i9*(T% s e hS )}9tJ kp N k }p k+1 ...p N 



hi 

h=l 
N 



= £>i...p fc _i[-40W+ [V k ,T a ] a p[fJ k (r) ■ A a± (r,fJ k (r)) - 

k=l 

-g] J d 3 aK ab (rf k (T),a; r)pf M) (r, a)}p k+1 ...p N + 

+ (terms containing factors N% ) = 0, (127) 

so that, with iV| = for every h, the equations for v ka {r) are 

dVk ^ - c abcV k c(T)[fJ k (r) ■ A b± (r, f) k {r)) - 

- gl J d 3 aK ab (Ur),a; r)pf M \r, a)] 4 0, (128) 



in accord with the mean value of the last of Eqs. ( |108| ) and with Ref . [ I3| except for the last 
term which is absent when the field is considered external. 

In analogy with the electromagnetic case (see Ref. |15||), we could try to find so- 

lutions to the field equations ( [L17| ) necessarily of the form A a ±_ (r, a) = A a± (r, a) + 

•Uo) 

(terms at least linear in the Qi a (T)' s) [A a _ L (T,a) satisfies Eq . (\1 1 7|) with Qi a (r) = 0, namely 



in absence of particle sources], to put them in the particle equations ( |118| ) and to take the 
mean value (< D >= J pDdp) of the resulting equations to get the real "classical" equa- 
tions. In the electromagnetic case, this procedure eliminates infinities from self-energies and 
causal pathologies (runaway solutions or preaccelerations) from the classical equations (see 
also Ref. JTJI for the rest-frame analysis of these problems) , which, instead, would be present 



by taking immediately the mean value of Eqs. flllTp and ( |118| ). However, it is not possible 
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to attack the former procedure in the non-Abelian case in absence of an analogue of the 
Lienard-Wiechert potential [but see in any case next Section for the N=2 case]. 
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VI. THE QUARK MODEL AND THE PSEUDOCLASSICAL ASYMPTOTIC 

FREEDOM FOR N=2 



In the nonrelativistic quark model, in which no SU(3) color Yang-Mills field appears, 
one assumes that the physical states are color singlets. Since Eq. ( |107|) gives the invariant 
mass of N colored relativistic scalar particles together with the SU(3) color Yang-Mills field, 
this is the starting point to try to extract a pseudoclassical basis for the missing relativistic 
quark model. 

A first step is to study what happens if we add 8 extra constraints implying the vanishing 
of the total color charge, so that only global color singlets are allowed for the particles+Yang- 
Mills field system: 

N 

Qa = Q ( a M) +J2Q™ = °- ( 129 ) 
i=l 

Then, we ask that these 8 conditions be fulfilled by the separate vanishing of the particle 
and field contributions to the color charge 

i 

gi yM) = o. (130) 

The first condition defines the relativistic scalar quark model: the particles by themselves 
are a color singlet independently from the SU(3) color field. For N=l, Eqs. (|130|) plus the 
constraint N = J2 a 9* a ®a ~ gives 9 conditions on 6 Grassmann variables: therefore a single 
pseudoclassical scalar quark cannot be a color singlet. For N=2, besides the constraints 
Ni = J2 a 6*Jia « and N 2 = J2 a 0* 2 J 2a ~ 0, one has 

Qla + Q2a = 0, (131) 

namely 10 conditions on 12 Grasmmann variables. 

The condition Q^ M ^ = replaces the Abelian condition A± (r, a) = ti± (r, a) — of 
absence of radiation [see Ref. J7|]. Since in the non- Abelian case we do not know how to 
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solve the equations of motion and since the superposition principle does not hold, we can 
only ask that there is no color flux on the surface at space infinity. This requirement also 
implies that the pseudoclassical solutions of the SU(3) Yang-Mills equations are restricted 
to those configurations which are color singlets like the glueball states at the quantum level. 

Due to the Gauss laws, the condition Q a {.T~) = can be imposed by choosing suit- 
able boundary conditions on the transverse SU(3) fields (see Eqs.(2-40) of Ref. (H): 
A a ±(r, o ! )^| ( ?|^ooO(|ap (2+e) ), 7r a± (T, o ! )^| (? |^ 00 0(|ap (2+£) ) with e > (for e — > one gets 
Qa 7^ 0). With these boundary conditions the requirement of color singlets for the whole 
theory becomes automatically the same requirement in field-independent quark models, if 
the equations of motion imply Q a YM \f) = 0. 

Eqs.(|l30D, (|IDD and ([109]) imply 



4-(Qia(r) + Q2a{T))=c acd Q ld (r){rj l {T) ■ A c ±(r, 7^ (t)) - ff 2 (r) ■ A c± (r,ff 2 (r)) - 
-g 2 s J ^[Kebifhir), a; r) - K^r), a; r)]pf M) (r, <?) + 
+9s[K c b(fJi(r),ff l (T);T) + K cb (ff 2 (r),f] 2 (T);T) - 

-K c6 (7f 1 (r) ) %(r);r)-K c6 (?f 2 (r),77 1 (r);r)]g 16 (r)=0. (132) 

In the gauge A(r) = g{r) = 0, the equations of motion for the field and the particles and 
the equations defining the rest-frame become [in Eq. (|133|) there is the same F-function of 
Eq-Ol 



Pl s {B){5 st [5 ad {^ + A) + c ad J c± (r, <})—]■ 

K uv (a, h t)D^\t, B)]^- }A d± (r, B) ± 
= ~ 9 2 s Pl s ^){Qia(rM(r)5 3 (a - Hfc(r)) - 
-r] s 2 (r)5 3 (a - ff 2 (r))} +c abd A' dL (r, a) J d 3 a' K bc (a,a ;r) 
W M \r, a') - Q lc (r)[5 3 (a' - ft(r)) - 5\a' - ff 2 (r))}} 
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-g 4 s P r AB)A s d± (r,B) J ' d 3 a,d 3 a 2 

[p^V^i) " <Ur)[*Vi " fh(r)) - *Vi - ff 2 (r))} 
[F a bcd{v ,Vi,<?2',t) + c auv c vde K eb (a ,&i;t )K uc (a ,a 2 ;r)] 
[p? M) (r,a 2 ) - Q lc (r)[5 3 (a 2 - ffi(r)) - 5\a 2 - fc(r))], 



d r 1)1 

- r [r] 1 m 1 



drl ' ^Ji-m(r) 

= T,QUr){KAr,Vi(r)) + fe(r) x M a (r, 771 (^H - 

a 
a 

J d 3 aK de (ff 1 (T), a; r)c enb A n± (r, a) ■ E b ±(r, a) + 
+^E{0i«(r)0i6(r)[Pa*A.(r,^(T))(ir c6 (^(r),^(r);r) - 

-X cb (?7i(r), 772 (r);r) ) + 7^l<?=?n ( K ab(v,f]i( T ); T ) ~ ^a&(<?, ^(t); r) )] 



9X a6 (77i(r),CT;r) (yM ) 
jp fa (r,(j)}, 



d?7i 



= " E0ia(r){4\(r,r/ 2 (r)) + fe(r) x 4(r, t^t))]'} + 

a 

+ E Oia(r)Pr(?7i)camA±(r, 7fc(r)) 

a 

y d 3 aK de (ff 2 (r), a; T)c enb A n± (r, a) ■ E b± (r, a) - 
-5 s 2 E{^i«( r )^( r )[ c ^d±(r,^2(r))(i^ c6 (77i(-r), ?7i (t);t) - 



a, 6 



9 

Kcbfair), 772 (r); r) ) + t^|<?=^ (^a&(^, »7i(r); r) - if a &(^, 7/2(7); r ) )] 



9a 

-Qio(t) / d 3 a{c adc A d± (T,rf 2 (T))K cb (rf 2 (T),a; r) + 



dK ab (f) 2 (T),a] r) 1 „ ( y M ) 



^7/5 

ryimi /U J = + 77 2 m 2 = - 
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-J2Q la (r)[A a± (r,fi 1 (r))-A a± (T,fj 2 (r))+g; 2 f d 3 a{dA s a± E s aX }(T,B) = 0. (135) 

Note that in Eqs.( |134|) the Coulomb interaction inside the kernel K does not contribute to 
the terms Qi a (r)Q lb (T) due to c adb Q la Q lb = S ab Q la Q lb = 0. 

The invariant mass of the system is [we choose r/i = 772 = +1, i.e. the quark and antiquark 
have positive energies and the antiquark is distinguished by the opposite color charge as a 
classical antiparticle moving forward in r (see Ref. |]20|-p2| and its bibliography)] 



+ Jm\ + (£ 2 (r) - Qia(r)i±(r, r/ 2 (r)) ) 2 + 

+ \y%, Xl, fa±](r) + dM9sK± + 9; 2 B al _](r, a). (136) 

In the case N=2, which should correspond at the quantum level to a meson configuration 
formed from a scalar quark and a scalar antiquark, plus glue [there are not sea-quarks, 
because there is no pair production at this pseudoclassical level], if N\ — > X^=o ^L^ia — 1 
and N 2 —>■ J2a=o^2a^2a — 2, the first condition in Eq. (|130| ) and the observation that the 
SU(3) Yang-Mills fields, solutions of the field equations ( |133| ), will depend on the Grassmann 
variables of the particles only through the color charges Qi a {T~), allow us to write the following 
developments [from now on we shall not write J2 a f° r repeated color indices] 

A*l(t, B) = A a± (r, B) + Q lu (r)A au± (T, B) + Q lu (r)Q lv (r)A auv± (T, a), 

7?a±(r, B) = vr^(r, B) + Qi«(t)7T 1 ± (t, B) + Qiu(t)Q 1v (t)ti { ^ v± (t, a). (137) 

As a consequence we have [we suppres the r-dependence and also the a-dependence when 
possible with the replacement Bi — ► i] 

2 



pf M M(B) + Q lu pZ Mm (B) + Q lu Q lv pZ^ 2) {B) + J> 

1=1 

r,(o) x ~Jfl) , M) ~M) *(«,) 

Cabc[A b± ■ 7T C± )((T) + Qiu[C abc {A b± ■ 7T CU± + A bu± ■ 7T c± )(ct) - 
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o- 



-S au (S 3 (a-ff 1 )-S 3 (a-ff 2 ))} + 

+ <3iu<3iuC a &c[A>± • + " k cv± + A buv± ■ 7r cJ J(cr), 

Pa{Pl)pb{P2) = C amn A^ ± (l) ■ 7f2(l)c 6rs l^(2) • 7?2(2) + 

+ Qiu[(c amn (Al(i) • #!la) + XLw • - 

- <W<s 3 (*i - ft) - ^(<?2 - m)) )c brs A r %) • #2(2) + 

+ c amn A^{l) • #2(l)(c fers (Xt(2) • #12(2) + Xl(2) • #2(2)) - 

" S bu (5 3 (a 2 ~ Vl) ~ ^{$2 - ff 2 )) } + 

+ QiuQiAc am nA m± (l) ■ <{(l)c 6rs (A r± (2) • < u 2(2) + 
+ XI (2) • #2±(2) + XL (2) • #2(2)) + 

~ (2) ~-J^) - (1) 

+ (1) + 

+ XLw • #2a)K s Xt(2) • #2(2) + 

+ ( c amn (A m± (l) • <2(1) + A mu± (l) ■ <|(1)) 

- -M^Vi - ft) - <* 3 (<?i - ft)) ) 

■ ( c brs (Xt(2) • #2x(2) + XI (2) • #2(2)) - 
-S av (5 3 (a 2 -ff 1 )-5 3 (a 2 -ff 2 )))} = 

= R(l,2)ab(<?l, °2\t) + R(l)ab(&l, &2] t) + R(2)ab(&l, <?25 7") + Rab(&l, &2',t), 

R(i,2)ab(vi,V2;r) = Qi a Qib[5 3 (ffi - ft) - 5 3 (ai - ff 2 )][5 3 (a 2 -771) - 5 3 (ai -572)], 

R(i)ab(&l, °2] r) = (-yQlu[$aJ 3 (Vl ~ Vi)c brs A r± (2) ■ 7r2(2) + 

+ 5 bu 5 3 (a 2 - ffi)c amn A m± (l) ■ 7f2(l) = Oiu(T)-R(i)at>u(tfi, ^2; r), 
# a &(<?i, ^2; r) = R$(#i, o 2 \ t) + Q lu R { ^ u {a u a 2 ; r) + <3iuOi«-R2L(^i> ^2; r) = 
= c am „Xl(l) ' #2(1)0™X°±(2) ' #2(2) + 

+ Qiu[(c amn (Al(l) ■ #H(1) + XL(1) ' #2(1)) )^X°±(2) ' #2(2) + 
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+ c am J^ ± (l) ■ ^i(l)(Q, rs (Xt(2) ■ 7^(2) + 1^(2) • #2(2)) ] + 

(2) + 

+ Xl±(2) • #2±(2) + XL (2) • #2(2)) + 

%(°) - (2) "-i 1 ) - (1) 

+ 

((4»±(i) (1) + A mu± (l) (1) + 

(l)) Cbrs A r± (2) ■ nll(2) + 

+ ( c amn (2^ ± (i) • #^(1) + ^(i) • #2(1)) ) 

• ( c 6rs (l2(2) • #2x(2) + Xlx(2) • #2(2)) ) ]• (138) 

In the last lines we separated the particle-particle -R(i,2)a&], the particle- field [R(i) a b] and 
field-field [R a b] terms. 
If we write 

Clf ^(ai, *2j r) = (if Vi, <? 2 ; r) + giu^'^ Vi, * 2 ; r) + 

_ 4(1) 4(2)-, 

+ " " r), (139) 

where (if \ff lt a 2 ; r) = c(a, - a 2 )(Pe^ d ° ^^"^ ) ab and and (^'^'^ 

are functions which could be evaluated by using the definition of Wilson path-ordering, then 
we can write the following decompositions of the interaction kernel K ab (ai, <j 2 ; t) 

Kab(ei, r) = K^(&!, a 2 ; r) + Q lu K^(^i, o 2 \ r) + Q lu Q lv K^ b ] uv (a l7 a 2 ; r), 



SabS 3 (a 3 - a l )5 3 (a A - a 2 ) 



Aix\a 3 - (T 4 | 



+ 



-Jp) 



_ Pfa - ff 2 )[A"(T, ff 3 ) • C^ 0)) Uff 3 , gi; r) + gi <-> ff 2 ) 

1 1 — » — » 1 



+ 



47r|a 3 - <t 4 | 

[l2( r > 3») • C (A<l0>) (^3, ^i; r)] au [l2( r > ^4) ' C^H^, B 2 \ r)] bu 



d 3 a 3d 3 a 4 
4n\a 3 - <t 4 ] 



Jo) 



An\a 3 - <t 4 | 



[5 3 (a 4 - a 2 )[A ± (r, a 3 ) ■ ' ^ \^a l -r) + 



•ai) 



+ £u!(T, <7 3 ) ■ C ( ^ D)) (^3, <?l! T) ] ab + (<Ti <-> <7 2 ) + 
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+ [A ± (r, as) • C {A - Vs, Bi\ r) } ae [A ± (r, <? 4 ) • ^ ^ >(* 4 , a 2 ; r) + 
+ K ± {t, «t 4 ) ■ C ( ^ 0)) (^4, <? 2 ; r)]* + 

+ [A ± (r, * s ) • CT- ' " } (<? 3 , ft; r) + A„ ± (r, ft) • ( {A - \a 3 , ft; r)] ae 

[A^\r,a 4 ).C (A -\a 4 ,a 2 ;r)] be }, 
~(2) / ~* ~* \ f d 3 a 3 d 3 a 4 

4°) -W(°) iW 1(2)n 4 1 ) 

■ [ 5 3 (ft - ft) [A ± (r, ft) • ^ * ^ } (ft, ft; r) + A u± (r, ft) ■ 

• 5 X ,A± ] (ft, ft; r) + A uv± (t, ft) ■ >(ft, ft; r)] ab + (ft <-> ft) + 

4°) ->, J(VK 4°) -^(°) J( j )l 

+ [A ± (r, ft) • >(ft, ft; r) } ae [A ± (r, ft) • & ^ } (ft, ft; r) + 
+ A u± (t, ft) ■ ^ ' ± '(ft, ft; t) + A uv± (t, ft) ■ } (ft, ft; r)] be + 

+ [A ± (r, ft) ■ (i A - >(ft, ft; r) + A u± (r, ft) ■ C (A ^(ft, ft; r)] ae 

4°) -iA (o) A {1) ) 4 1 ) -}!(»), 

[A ± (r, ft) • Q x ' x (ft, ft; r) + A„ ± (r, ft) • C (Ai '(ft, ft; r)] be + 

+ [A ± (r, ft) ■ <^ ' A± ^ Vs, ft; r) + A u± (r, ft) ■ ^ ^ \a 3 , ft; r) + 

+ A~2 ± (r, ft) ■ C^Hft, ft; T)] ae [l { ° ] (t, ft) • £*(ft, ft; r) ] be ], (140) 

and of the potential V in particle-particle [Vpp], particle- field [V(i)pp] and field-field [Vpp] 
contributions 

1 - -4o) 2 -4°) -( ) 41) 

-V[ffi, A a± , 7f a± ](r) = Vpp[^, 77^ A a J(r) + ^ V^pp^, A a± , 7r a± , A ouJ _](r) + 

z i=i 
+ V FF [A ± ,7?r](r), 



Vpp = 9 2 s J d^aid^^Y^Rii^abKa^ax.a^T) = 

a,b 

= g 2 s J d 3 a l( i 3 a 2 £ OiaOitl^^i " " 5 3 (^i " ^)] 

a, ft 

[ 5 3 ((?2 _ ^) _ 5 3 ((?2 _ rf2 )]^) (c?1) a,; r), 
V(i)p F = 9 2 s J d 3 aid 3 a 2 Y^[R(t)abK ah ](a u a 2 ;T) = 

a,b 

= 9s d 3 a x d 3 a 2 ^ Qi u R(i)abu(^i, ft] T)[K { ° b ] + Q lv K^ v }(a u a 2 ; r), 
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V FF = g 2 s J d z a x d z a 2 p^ M \r, Bx)K ah {a x , a 2 ; r)/f A/) (r, 30 = 

a,b 

+ Q lu Qiv{R^Kfl v + R { lK^ + R { l v K<$)]{a u a 2 ; r). (141) 



The explicit form of the particle-particle potential is 
Vpp = \glY,QiaQib J d 3 a, J d 3 a 2 [5 3 (a 1 -ff 1 (r))5 3 (a 2 -ff 1 (r)) + 

a,b 

+ 5 3 (°i ~ m(r))S 3 (a 2 - ff 2 (r)) - 2S 3 (a 1 - m(r))5 3 (a 2 - ff 2 (r))} 

• -1 [d 3 * 3 [d 3 ad 6 ^ Z M^^ + 
An J J | 03 — 04 I 

'JO) ,j(o)s 

+ c au „ £ 3 (0 3 - 01) (04 - ^2) • A v± (t, ^4)Cb X (^2, ^4; r) 

4vr I 03 — 04 1 1 04 — 2 | 3 

-JO) fA {o) ) 
+ ^ 3 (04 - O2) (^3 - gl) • A^kfc g 3)Cfc X (gl, T) , + C arl ,C fes f 



^3 - <?4 || 0*3 - o x | 3 ' (4vr) 2 

-(°) (A (o) ) -(°) (A (o) ) 

(04 - 2 ) ■ A„ ± (r, 4 )Cr M x (02, 04; r)^ - 01) • A t± (r, 3 )Cs M x (di, <?3; t) 
I 03 - 01 | 3 | 03 - 04 1 1 04 - 2 | 3 



= ^ 2 £<2ia<3i& / d 3 ax J dV 2 [5 3 (0! -fh{r))8 z {a 2 ~ Vi(r)) + 

a,b 

+ 5 3 (ax - ff 2 {r))5 3 (a 2 - rf 2 (r)) - 25 3 (0i - m(r))5 3 (a 2 - ff 2 (r))} ■ 

-JO) (A <o) ) 

— [ d 3 a 3 [ d 3 ^ { Cauv f 3 ^ 3 - ^ ^ 4 ~ ^ ' Av± ^ B ^ ub ^ ^ T ^ 1 
An J J An | 03 — 04 1 1 04 — 2 | 3 

-_(0) r A (o)s 

+ (04 ~ ^2X03 - gl) • A V± (T, 3 )C6 X (glj ^3; T) + C arv C bst _ 
I 3 - 4 || 3 - 01 | 3 (47r) 2 

'JO) -40) ( ^(°)) 

(04 - 2 ) ■ A v± (t, a^Cru 1 - (02,O4;t)(0 3 - 01) ■ A t± (r, 3 )C^ X (0i,0 3 ;r) 2 
I 3 - 01 | 3 | 03 - 04 || 04 - 2 | 3 

The first term, with the instantaneous Coulomb interaction, does not contribute, because 
for every N we have J2 a QiaQia = 0, due to Eq. (|23|) . 

With regard to the other two terms, let us remark that in the limit 171 = 172 [and thus 
also in the points of maximal divergence, i.e. B\ = o 2 = for the first two terms and 
0i = 02 = = 0' for the third one] the potential vanishes due to the multiplicative term 
containing the Dirac functions. 
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With the substitutions a = <j\ + £ and a 1 = <?2 + the particle-particle potential can 
be put in the form 

d^ 



V PP = 9s H QlaQlb{ 



T>, 1 (4>r) 2 1 { I 2 



i=l 



cS%(r),*(r) + gr)- 



ell^(r)-^(r) + e 



\i\Ux{r)-^{r)-i 
1 /• rf 3 ^ 3 f 

~r Cart; Q>si 



C^(^(r)^ 2 (r)+e;r)] + 



3 



J (4tt) 3 I i | 3 | £' 
£ tMirMr ) + £)£> ■ Afl (r, ff^r) + £') 
i=i 



- 2 



£ ■ a { :1(t, m(r) + £)£> ■ a[°2(t, Ur) + £') 



c£°\fk(T),Vi(T) +£;r)C^ ) \ff i (r),v i (r) +£';t) + 



\fh(T)-7k(T)+Z-? 



• ciu l) \ffi(T),ff 1 (r) +£;r)C ( s ?\ff 2 (T),ff 2 (r) + £>)] = 
_ 9s[Vi -V2;Vi;A ± ](r) 

- ^2; A ± ](r) 0. (143) 

This is the pseudoclassical statement of asymptotic freedom for N=2. 
By using the developments ( P-37| ), the equation of motion for the color charge Q\ a {r) = 
-Q2.{t) becomes 

d ~ ■ 

-r-Qia(r) = CacdQuiT^ff^T) ■ A^faffifr) - 

-9 2 s J <P<rK*{jh.{T), d\ T) Chuv A u ] ± (T, a) ■ tt2(t, a)] + 

+c ac dQid(r)Qir(r)[ff 1 (T) ■ A cr± )( r ^i( r ) + 

+9 2 s { ft (r) , ft (t) ; r) - JT<? (m (r) , fc(r) ; r) } - 

-9 2 s J d z a{K { ^{f] 1 (T) : a; r) Cbuv (A u ] ± ■ t$ ± + aJ ± ■ 7r2)(r, S) + 
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+ K i2(Vi(r), a; T)c buv A u l(T, a) ■ |uec7r^(r, a)}, 



while the particle equations become 

d r){(r) o 

[mi- 



= Y,QUr){E { :l r {rMr)) + fo(r) x B°\r, ft(r))n - 



-<?7 dVfc*^* (r,^(r))^(^(r),?;r) + 



+ ^ )( t ( ; ) '" ;r) ]pr )(o) (r,^)}+ 

(777! 



+ EOiaWQi tt (r){^ 1 i r (r,r/ 1 (r)) + [ft(r) x (r, ft(r))f - 

au 

~_» ^(1) -41) ^(o) „, . , . 

c enfe (A n± • + • £ 6± )(r, a) + (C(r, m(r))^(77i(r), a; r) + 
+i« s ± (r, r/ 1 (r))^ ) (r/ 1 (r), a; r)) Cen6 (ll°i • l^)(r, a)} + 



+S S 2 [c a ,c4? (r, ft (r) ) (tf £> (ft (r) , ft (r) ; r) - £> (ft (r) , ft (r) ; r) ) + 
+ ft( T ) ;T ) _ X (o)(^ft (r);r))] _ 

+ (c arfc ( (r, ft(r))^(ft(r), a; r) + a£1(t, ft (r))l^(ft (r), r)) + 

, di^(ft(r),^^( yM)(o) ^ j 



% 2 - ^ (r) 



= - £0ia(r){^? (r,ft(r)) + [ft(r) x B? (r, ft(r))H - 



-Pr(ft)c amrf i^ s (r,ft(r)) I ^^(^M.ffirJce^Cr,?) -^(t,*) 
-^ 2 / rfV[c adc iif (r,ft(r))^i o) (ft(r) !( ?;r) + 
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i dK^(ff 2 (T),a;r) { ym)(o) ( ^, 

an 



~4o) ^(1) -41) -_(o) „, . . . 

c enfe (A n± • + A nu± • £ b J(r, a) + (^i s (r, ^^^(^r), a; r) + 

fk(T))K£\ih(T), a; T))c enb {Zl ■ E b °l)(r, a)} + 
+g 2 s [ca d cA { d T(T, fj 2 (r))(K^(ff 2 (r) : m(r); r) - K^(ff 2 (r), ff 2 (r); r)) + 
+ ^\a=r} 2 (K^(a,m(r);r) - K^(a : ff 2 (r)-r))} - 
-9lj d 8 (7[p B ^ r (^*(r))^ ) (^(r),?;r)p£: M 5«(r,?) + 
+ (c arfc ( iif (r, ff 2 (r))K^ff 2 (r), a; r) + i£l(r, m(r))K^mr), a; r)) + 
+^%|^]pr )(o) (r,^)]}, (145) 

and the equations defining the rest-frame are 



mi 1 + m 2 



-4°) -4°) - -41) -41) 

-J2Q^(r)[A a± (r,Vi(r)) ~ A al _{r,fj 2 {r)) + Qi„(r)(^ ou± (r,7fi(r)) - A au± (r,ff 2 (r)))} + 

a 

+Qiu(r)Q lv (r)(dAif E^ ± + dA^E^l + dA^ ± E^ s )}(r, a) ± 0. (146) 

The invariant mass (the relative Hamiltonian) of the pseudoclassical scalar quark model 
takes the form 

/ " -40) -41) T 2 / I -40) -41) T 2 

H rel = ]JMf[m 1 ,fl 1 ,K 1 ,A a± ,A ± ](r) +Ki(r) + V M|[m 2 , ff 2 , k 2 , A a± , A ± ](r) + k 2 (t) + 

-40) 2 -40) „ (o) -41) 

+ Vp P [?7i - 772; 4J( r ) + X! V(i)PF[(#; 4u_; ^±5 4l 10") + 

i=i 

+ V> F [i fc) , \r) + \jd*a Y:tfX ± + ^ 2 ^](r, a) (147) 



with 
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-Jo) ~Jl) 
Ml = m\ + R m 2^(r) • Y.QUr)[A aX + £ Qiu(r)A au± }(r, #(t)) + 

a u 
•JO) -Jo) 

+ J2QUr)Q lb (r)[A a± -A b± }(r,ff l (r)). (148) 

a,b 

The first three terms of H re i, after a suitable average over the field degrees of freedom, 
should give the effective rest-frame Hamiltonian for the pseudoclassical relativistic scalar 
quark model: while are the quark current masses, the suitable average of M, should give 
the quark constituent masses. The terms Vu\pp describe quark-field interactions. The last 
two terms should describe the pseudoclassical glueball degrees of freedom. 

All the previous results can be reexpressed in terms of rescaled Yang-Mills fields by 
putting A aJ _ = g s A aJ _, 7r a± = gJ l TT a ^ B a = g~ l B a - 

Finally, the Berezin-Marinov distribution function is 

P = P1P2, 

Pi = -2 < Q la {r) > Qia(r)Ni + ~iVf « 0, (149) 

6 

and, by taking the mean value of equations with it, the "classical equations" are [see Refs. 
|T5| , |23|p4]| for the case of an external Yang-Mills field] 

d - "J°) 

< QUr) > -c abc < Q lc (r) > [^(r) • A b± (r, ff^r)) - 

-9 2 S J dV^( 7 f 1 (r)^;r)pr M)( ° ) (r,a ; )] = (150) 



dT yfi - m(r) 

= E< QUr) > {Eif(r,ffi(r)) + fe(r) x B°\r, ff 1 (r))Y} - 



P r ± s (ffi)c amd A^'l(T, ffi{r)) J d^aK^ifftir), a; T)c enb A n± (T, a) ■ E b± (r, a) 
g 2 s J dMca d cA { ;l r (r,Ur))K^mr),a-,T) + 

}pr )i0) (T,a)}, 



dr]\ 
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d -[m 2 ^ 



<,Tl ""y/i-v 2 (r) 



= ~ E < QUr) > {E { ;l r (r,Ur)) + Mr) X B^\r, fJ 2 (r))Y} 



-Pr^CarnAt^Ur)) J d 3 aK { d °J (fj 2 (t) , a; r)c enb A°l (r, 3) • I SJ(r, a) 
-£ 2 / dV[c odc i(f (r^ 2 (r))^?(%(T),a;r) + 



^ )(?72(r) '^]pr M)(o) (r,-)}, (151) 



Vl\T) , V2\. T ) . -2 f .3 /^(o)iA( ) lv 

™i / . 2 = + ra 2 ^== + g s da{dA)J_E^l){T,<T)- 
Vl-^(r) Vl-%(r) J 

-E < Oia(r) > [^£W(r)) - 1^ (r, fc(r))] + 

a 

E < QUr) > ! dMdAifE^ + SAgktf* )(r, a) = 0. (152) 

a, it J 

The particle "classical equations" are expected to have the same causal pathologies (run- 
away solutions or preaccelerations) as the Abraham-Lorentz- Dirac ones in the electromag- 
netic case. In that case, see Refs. [0,0 , one can follow a different procedure: i) solve 



the field equations (the solutions are incoming free linear waves plus the Lienard-Wiechert 
potential if retarded solutions are selected); ii) put the solution in the particle equations; 
iii) take the mean value of the new equations. This procedure gives different "classical 
equations" without pathologies, because one has Q\ = for the Grassmann-valued electric 
charges: in this way, at the pseudoclassical level, one regularizes the Coulomb self-energies 
0, has no radiation produced by the single Lienard-Wiechert potentials but has radiation 
from the superposition of them [JTJJ] (terms in QiQj with i ^ j). 

In the non-Abelian case, we do not know solutions of Eqs.( |133| ) and the superposition 



principle does not hold, so that we do not know how to apply the second procedure. A first 
needed step would be to find the Green function of the operator on the left side of Eq.( |133| ). 
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VII. DISCUSSION AND CONCLUSIONS. 



We have obtained the description of the isolated system of N scalar particles with 
Grassmann-valued color charges plus the color SU(3) Yang-Mills field on spacelike hypersur- 
faces and then its canonical reduction to the Wigner hyperplane with only physical degrees 
of freedom (generalized Coulomb gauge in the Wigner- covariant rest-frame instant form of 
dynamics) following the same steps of Refs. |J7,n| for the case of N scalar particles with 
Grassmann-valued electric charges plus the electromagnetic field. 

In particular, we obtained the reduced Hamilton and Euler-Lagrange equations for 
the reduced transverse Yang-Mills field and for the particles in this rest-frame generalized 
Coulomb gauge, generalizing those of the electromagnetic case. Four kinds of interactions 
between two color charge densities (associated either with a particle or with the SU(3) 
field) were identified in the resulting potential extracted from the gauge part of the Yang- 
mills field: one of them is a Coulomb interaction, while the others involve Wilson lines 
[C(<?i> &2', t) = Pe^ da Aa - L ^ T ' a ^ Ts j a i on g straightlines between the two (simultaneous in the 
rest frame) points where the two color densities are located. In the electromagnetic case 
|7Jl^], where the field is not charged, there is only the interparticle Coulomb potential. Due 
to the high nonlinearity of the field equations, it is extremely difficult to try to define a non- 
Abelian analogue of the Lienard-Wiechert potentials of the Abelian case MM. Moreover, we 



have found the Berezin-Marinov distribution function for the Grassmann sector |T5||T6] and, 
then, the "classical equations of motion" as the mean value of the previous equations. In 
absence of a Lienard-Wiechert potential, solution of the field equations, to be put in the 
particle equations of motion (which, otherwise, are going to have causal pathologies like in 
the electromagnetic Abraham-Lorentz-Dirac equations), we cannot take the mean value of 
these new equations and check that the resulting equations do not have causal pothologies 



as it happens in the electromagnetic case [14 



We then studied the N=2 case, which is the pseudoclassical basis of the relativistic 
scalar-quark model but with the reduced transverse color field present, which describes the 
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pseudoclassical glueball degrees of freedom. With suitable constraints on the Grassmann 
variables we obtain the description of an isolated system corresponding to a meson: a quark- 
antiquark pair plus the transverse SU(3) field. The reduced Hamiltonian of the system is 
its invariant mass expressed in the intrinsic rest frame. It contains the relativistic kinetic 
terms of the two quarks minimally coupled to the transverse SU(3) field, an interparticle 
field- dependent (but quark-mass- independent) potential, a single particle - field potential 
and the kinetic terms plus a self-interaction potential for the transverse SU(3) field. If 
one would know the non-Abelian Lienard-Wiechert potentials of the quarks (if this concept 
makes sense due to the nonlinear selfcoupling of the color field, i.e. due to the glueball degrees 
of freedom), one could obtain a field- independent expression for the invariant mass of the 
meson and therefore an explicit formulation of a field-independent relativistic scalar-quark 
model starting from the pseudoclassical Lagrangian for QCD with scalar quarks. This would 
fill the gap, at least at the pseudoclassical level, between QCD and quark models [see Ref. [13| 
for the nonrelativistic one]. Moreover, a field-independent (but quark- velocity-dependent) 
interparticle potential would appear, whose static [m, — > oo =>- 77$ (t) — > 0] part should have 
connections with the static potential in QCD P5fl . Note that at the pseudoclassical level 
there is no notion of gluons (or gluon exchange): only SU(3) fields are present and only a 
Lienard-Wiechert potential (if it exists in some sense) would create a bridge towards the 
Wilson loop expectation value form of the static potential between heavy quarks [in the 
quenched approximation in which quark loops from pair production in vacuum (absent in 
the pseudoclassical theory) are neglected (sea-quarks of infinite mass) and glueballs (closed 
color loops) are unambiguously defined] re-expressed in terms of perturbative QCD. 

When the study of Dirac fields and spinning particles on spacelike hypersurfaces will be 
finished [|12], one will get analogous results for spinning particles (with Grassmann- valued 
spin), namely one will introduce the quark spin structure in the pseudoclassical relativistic 
quark model. 

In the N=2 (meson) case, we have explored the implications of the imposition of the 
condition that the isolated system (quarks+transverse SU(3) field) is a color singlet. The 
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condition Q a = Qa ( T ) + Qi (r) + Q2o,(t) = is imposed by hand and, moreover, it is 
asked to be fulfilled by asking separately the two conditions Oa ( r ) = [ no c °l° r nux of 
the transverse SU(3) field at space infinity; it replaces the Abelian condition of no radiation 
field, i.e. A±(r, a) = [0,0]] and Qi (r) + Q2a(T~) = [it is the color singlet condition of 
field-independent quark models, which work very well phenomenologically |13|]. Now, the 
condition Qair) = can be imposed by choosing suitable Hamiltonian boundary conditions 
at fixed r on the transverse SU(3) field. Therefore, if confinement exists, the condition 
Q^I M \ T ) = Qia( T ) +Q2a(T) = should emerge from the solution of the Hamilton equations 
[at least approximately seen the phenomenological soundness of the quark model]. Moreover, 
the condition Qi a (r) + Q2a{i~) = 0, together with the other constraints on Grassmann 
variables, implies the disappearance of the Coulomb term in the interparticle potential [this 
is a statement stronger than the regularization of the Coulomb self-energies in the Abelian 
case R|l4l ] and that this potential tends to zero when the two quarks tend to the same 
spatial position in the rest frame [|^i(t) — ^(t)! 0]. This is the pseudo classical statement 
of asymptotic freedom: there is an antiscreening even stronger than in QCD, because the 
reducing (screening) effect of pair production is here absent. This kind of asymptotic freedom 
is an algebraic consequence of the request of color singlets in the quark model oriented form 
Qia(r) + Q 2 a(r) = 0. 

It would be interesting to study in this way \Qi a (j) + Q2a(j) + Q^aij) = 0] the N=3 
(baryon) case of 3 quarks or 3 antiquarks. What happens when either only two quarks or 
all three quarks tend to the same spatial position? 

Coming back to the N=2 case, the main unsolved problem, connected with the color 
singlet requirement, is confinement. Even if the boundary conditions for the transverse 
SU(3) fields are chosen so to imply Q a {r) = 0, we do not know how to do either analytical 
or numerical calculations to check whether the interparticle potential implies confinement 
and, if yes, whether confinement implies Qi a (r) + Q2a{ T ) — 0, not to speak of the possible 
glueball degrees of freedom. The obstruction is the lack of control on the Wilson line 
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C(<7i> c?2; t) = Pe^s dfT Aa± < - T ' "- ) ^ s ; given the function space for the transverse SU(3) potential 
^4a±( r , 0% which is the behaviour in ax and cr 2 of ^2! r ) as an element of the group SU(3) 
in the adjoint representation? Does it belong to the same function space as the transverse 
SU(3) potential? How to simulate it on a lattice having eliminated all the gauge degrees of 
freedom in favour of a transverse potential? Which is the pseudoclassical analogue of the 
quantum Wilson criterion of confinement p6j|27l in this Hamiltonian generalized Coulomb 
gauge? Moreover, an aspect of the confinement problem which has to be understood at the 
pseudoclassical level is the role of the center Z3 of SU(3), namely the zero triality condition 
[fermions know SU(3), but the Yang-Mills field feels only the not simply connected group 
SU(3)/Z 3 relevant in the approaches of Refs. p8| . Following Ref. [P9[] , this condition is 
probably hidden in the fact that the Wilson line C(o\, &2', t) must take values in SU(3)/Z 3 
and not in SU(3). 

Let us also note that the results of Ref. fl[ on Yang-Mills fields plus Grassmann-valued 
Dirac fields, once the description of Dirac fields in the rest-frame instant form will be termi- 
nated Hl2j ) suggest that the potential V will remain unchanged except for the replacement 
Ei Pia(r, a) i-> ^(r, a)T a ^p(r, a). 

Another yet unsolved problem (also in the electromagnetic case) is how to eliminate the 
3 constraints 7~C p (t) ~ defining the intrinsic rest frame. This requires the introduction 
of 3 gauge-fixings identifying the Wigner 3-vector describing the intrinsic 3-center of mass 
on the Wigner hyperplane. However, till now these gauge-fixings are known only in the 
case of an isolated system containing only particles. When the center of mass canonical 
decomposition of linear classical field theories will be available (see Ref. |3D| for the Klein- 
Gordon field), its reformulation on spacelike hypersurfaces will allow the determination of 
these gauge-fixings also when fields are present and a Hamiltonian description with only 
Wigner-covariant relative variables with an explicit control on the action-reaction balance 
between fields and particles or between two types of fields. 



As said in Ref. |fL4| , the quantization of this relativistic scalar-quark model has to over- 
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come two problems. On the particle side, the complication is the quantization of the square 
roots associated with the relativistic kinetic energy terms. On the field side, the obstacle 
is the absence (notwithstanding there is no no-go theorem) of a complete regularization 



and renormalization procedure of electrodynamics in the Coulomb gauge: see Refs. pl|j32 
for the existing results for QED. However, as shown in Refs. 0,f{]|§ , the rest-frame instant 
form of dynamics automatically gives a physical ultraviolet cutoff: it is the M0ller radius 



p = \ / —W 2 c/P 2 = \S\cfVP 2 (W 2 = —P 2 S 2 is the Pauli-Lubanski Casimir), namely the 
classical intrinsic radius of the worldtube, around the covariant noncanonical Fokker-Price 
center of inertia, inside which the noncovariance of the canonical center of mass x M is con- 
centrated. At the quantum level p becomes the Compton wavelength of the isolated system 



multiplied its spin eigenvalue ys(s + l) , p i— > p = J s(s + 1)%/M = ys(s + 1)Xm with 
M = y/P 2 the invariant mass and A m — h/M its Compton wavelength. 

Let us remark that in the electromagnetic case all the dressings with Coulomb clouds [of 



the scalar particles and of charged Klein-Gordon fields in Ref. [[14] and of Grassmann- valued 



Dirac fields in Ref. [TJ] are done with the Dirac phase r] em = —jrd- A ||. The same phase is 



used in Ref. |32|] to dress fermions in QED. In these papers there is a definition of dressing 
of fermion fields in the non-Abelian quantum case, which is quite similar to the one of Ref. 
HI] (used in this paper) even if implemented only perturbatively. Essentially, one looks for a 
matrix h G SU(3) such that under a gauge transformation U one has h *—>■ h u = U~ l h; then 
one has ip = hippuYS and A a T a = hA a puysT a hr x — dh br x with i/jphys and A a pHYS gauge 
invariant. Comparison with Eq . (|73|) and with Ref. |l] shows that at the classical level one 

Also in Ref. |29|] the solution of the quantum Gauss law constraint on Schroedinger 
functional *&[A] in the case of two static particles of opposite charges, is able to reproduce 
the Coulomb potential and the Coulomb self-energy with the same mechanism as in Refs. 
p|JT4]] [namely with the Abelian form of Eq.(|93|)] only if iff [A] = e iVem &{A] with gauge 
invariant and not with ^f'[A] = e dx A< ~ x ' x ^<f?'[A] with a phase factor resembling the Wilson 
loop operator [one has <&[A] = e dx A ^ x ' X ^<S>'[A], namely the Wilson line operator has been 
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broken in the gauge part plus the gauge invariant part using A = dr] em + A±]. However, it is 
difficult to see a connection between the phase of the Schroedinger functional *&[A] proposed 
in Ref. as a solution of the non-Abelian Gauss laws and our potential V, which is a 
consequence of Eqs. fl73|), flU, © [now the Wilson line e*C d£ -^ x "^ Ta cannot 

be broken in a gauge part and in a gauge invariant part due to Eq.([73|); in Ref. |29| the 
gauge-dependent part is a product U(xi)W (x ) of two gauge transformations]. 

Let us finish with some heuristic considerations about the M0ller radius. In QCD, due 
to asymptotic freedom and to the renormalization group equations, the strong coupling 
constant a s = g 2 /^ is replaced by the effective running coupling constant (see for instance 
REf. f3|) a s (Q 2 ) = 12tt/(33 - 2N F ) Irij^- at high Q 2 (N F is the number of flavors, 
giving the screening contribution of the fermions to the vacuum polarization). Therefore, 
the adimensional coupling constant a s may be replaced with the fundamental QCD scale 
(h — c — 1) Aqcd ~ 300Met> ~ 10 _13 cm = lfm, which is now usually replaced by a s (m 2 z ) « 
0.116 |H| [dimensional transmutation, connected with the breaking of scale invariance at 
high energies and with the scale anomaly; the physical mechanism for generating the scale at 
low energies is unclear (a candidate is the chiral symmetry breaking phase transition of QCD 
which generates a constituent mass of order 300 Mev for the light quarks)]. This implies 
|T3| that in the nonrelativistic quark model with confinement, one may choose (among 
the many possible phenomeno logical potentials) the simple potential V(r) = — |H^u _|_ Kr 
with the short distance behaviour a s (r) = 127r/(33 — 2Np) ln^%, i.e. for r < r Q = A( ^ CD ~ 
10~ 13 cm = lfm [in Ref. [25| it is shown that QCD perturbative results cannot be trusted for 
r < 0.07Aq^ d = 0.07/m for A QCD w 210Mev., i.e. a s (m 2 z ) = 0.118]. One can consider r Q as 
an effective radius of confinement for quarks and glueballs [the proton Compton wavelength 
is A p = h/rripC ~ 10 _13 cm = lfm < r Q ). For r w r G the expressions of a s (Q 2 ) and a s (r) 
break down due to confinement, which is described by the linear term in V(r) [k, « 0.2Gev 2 
is the "string tension", which turns out to be determined numerically as a function of Aqcd 
in lattice gauge theory |27j] at least for heavy quarks, for which a string-like (chromoelectric 
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flux tube) structure emerges; in the theoretical approach based on the analogy with type 
II superconductors (se the reviews in Refs. |34]Jrj|), where the vacuum is a color-dielectric 
medium and a qq state is a confined color-electric flux tube (anti-Meissner effect), k is 
determined by the gluon condensate 24 < J2a F a ^uF^ u > (as confirmed in strong coupling 
expansion of lattice QCD ||26||), which is present besides the < qq > condensate responsible 
for chiral symmetry breaking]. The crucial point for the pseudoclassical relativistic quark 
model would be to see whether Eq. ( |143| ) implies Vpp— ^-^i-kx, k \ffi — fj 2 \ + ■■■■ 



Let us also note that in the MIT bag model (see Ref. |35| for a review), the bag constant 
B is connected with the gluon condensate and the length of the chromoelectric flux tube for 
heavy quarks and the string tension k are determined by B and a s . 

The M0ller radius p = \S\/M is going to play the role of a ultraviolet cutoff p = 



s(s + 1)Xm [Xm is the Compton wavelength of the isolated system with invariant mass 
M = y/P 2 = H re i] at the quantum level (like the lattice spacing in lattice QCD). Since p 
describes a nontestable classical short distance region [impossibility of frame-independent 
determination of the location of the relativistic canonical center of mass (also named Pryce 
center of mass and having the same covariance of the Newton- Wigner position operator P5[]); 
its connection with the Mach's principle according to which only relative motions are mea- 
surable], it sounds reasonable that for a confined system of effective radius r Q = 1/Aqcd one 
has p ~ r Q . However, this is not correct because it implies a mass-spin relation \S\ ~ r Q M, 
while the phenomenological Regge trajectories are \S\ = a' s M 2 + a [a' s = IGev ~ 2 ], implying, 
at least for heavy quarks, an effective string theory inside QCD. Now, in string theory |37] 
the relevant dimensional quantity is the tension T s = l/2ira' s (the energy per unit length), 



which, at the quantum level, determines a minimal length L s = yh/T s = <J 27iha' s h = 1 ^2iia' s . 
For a classical string one has 15*1 < a' s M 2 , so that its M0ller radius is p < a' s M: at the quan- 



tum level one has p = \/s(s + 1)Xm < a' s M = L 2 S /\ M . Therefore, if the QCD string has 



L s = y2na' s ~ 10 13 cm = lfm < r Q = 1/ Kq C d ~ ^-f m , one gets that the M0ller radius of 
a confined system must be of the order p < r 2 M = M/ Aq CD [p < t^/Xm]- This corresponds 
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to a QCD string with 1txoi s <r 2 Q = Aq^ d : see Ref. [|38| for a long glue string giving rise to 
an effective Nambu-Goto string with this oc s and more generally see Ref. fl39f| . 

This effective QCD string theory (whose final formulation has still to be found) must 
not be confused with string cosmology [pTjMl], in which, at the quantum level, the string 



tension T cs = l/2irol cs = L 2 cs /Ti gives rise to a minimal length L cs h = ^2iia' cs > Lp [Lp = 
1.6 1CT 33 cm is the Planck length] and is determined by the vacuum expectation value of 
the background metric of the vacuum (if the ground state is flat Minkowski spacetime), 
while the grand unified coupling constant oiqjjt (replacing a s of QCD) is determined by the 
vacuum expectation value of the background dilaton field. This minimal length L cs > Lp 
(suppressing the gravitational corrections) could be a lower bound for the M0ller radius of 
an asymptotically flat universe, built with the Poincare Casimirs of the asymptotic ADM 
Poincare charges (see Ref. |4j] for the canonical reduction of tetrad gravity). The upper 
bound on p (namely a physical infrared cutoff) could be the Hubble distance cH~ x « 10 28 cm 
considered as an effective radius of the universe. Therefore, it seems reasonable that our 
physical ultraviolet cutoff p is meaningful in the range Lp < L cs < p < cH~ l . 

Let us remember that p is also a remnant in flat Minkowski spacetime of the energy 
conditions of general relativity [J7|: since the M0ller noncanonical, noncovariant center of 
energy has its noncovariance localized inside the same worldtube with radius p (it was 
discovered in this way) [Q, it turns out that an extended relativistic system with the 
material radius smaller of its intrinsic radius p has: i) the peripheral rotation velocity 
can exceed the velocity of light; ii) its classical energy density cannot be positive definite 
everywhere in every frame. 

Moreover, the extended Heisenberg relations of string theory f37j, i.e. Ax = + ^ = 

-t + T* 2 ( see Ref - for the 

meaning of Ap in string theory) implying the lower bound 



Ax > L cs = \jTi/T CS} have a counterpart in the quantization of the M0ller radius [Q: if we 
ask that, also at the quantum level, one cannot test the inside of the worldtube, we must ask 
Ax > p which is the lower bound implied by the modified uncertainty relation Ax = + 
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%3r. This would imply that the center-of-mass canonical noncovariant (Pryce) 3-coordinate 
z = y/~P*(x — -pzx°) cannot become a self-adjoint operator. See Hegerfeldt's theorems 
|43fl , his interpretation pointing at the impossibility of a good localization of relativistic 



particles (experimentally one determines only a worldtube in spacetime emerging from the 



interaction region) and also the comments of Ref. [44] against this interpretation. Since 
the eigenf unctions of the canonical center-of-mass operator are playing the role of the wave 
function of the universe, one could also say that the center-of-mass variable has not to be 
quantized, because it lies on the classical macroscopic side of Copenhagen's interpretation 
and, moreover, because, in the spirit of Mach's principle that only relative motions can be 
observed, no one can observe it. On the other hand, if one rejects the canonical noncovariant 
center of mass in favor of the covariant noncanonical Fokker-Pryce center of inertia P5|,|7|j4^[] 
Y^, {Y^,Y U } ^ 0, one could invoke the philosophy of quantum groups to quantize Y^ to 
get some kind of quantum plane for the center-of-mass description. 
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